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ABSTRACT 

We  study  Milky  Way  kinematics  using  a  sample  of  18.8  million  main-sequence  stars  with  r  <  20  and  proper- 
motion  measurements  derived  from  Sloan  Digital  Sky  Survey  (SDSS)  and  POSS  astrometry,  including  ~  170,000 
stars  with  radial-velocity  measurements  from  the  SDSS  spectroscopic  survey.  Distances  to  stars  are  determined 
using  a  photometric-parallax  relation,  covering  a  distance  range  from  ~  100  pc  to  10  kpc  over  a  quarter  of  the  sky 
at  high  Galactic  latitudes  (\b\  >  20°).  We  find  that  in  the  region  defined  by  1  kpc  <  Z  <  5  kpc  and  3  kpc  <  R  < 
13  kpc,  the  rotational  velocity  for  disk  stars  smoothly  decreases,  and  all  three  components  of  the  velocity  dispersion 
increase,  with  distance  from  the  Galactic  plane.  In  contrast,  the  velocity  ellipsoid  for  halo  stars  is  aligned  with  a 
spherical  coordinate  system  and  appears  to  be  spatially  invariant  within  the  probed  volume.  The  velocity  distribution 
of  nearby  (Z  <  1  kpc)  K/M  stars  is  complex,  and  cannot  be  described  by  a  standard  Schwarzschild  ellipsoid.  For 
stars  in  a  distance-limited  subsample  of  stars  (<100  pc),  we  detect  a  multi-modal  velocity  distribution  consistent 
with  that  seen  by  HIPPARCOS.  This  strong  non-Gaussianity  significantly  affects  the  measurements  of  the  velocity- 
ellipsoid  tilt  and  vertex  deviation  when  using  the  Schwarzschild  approximation.  We  develop  and  test  a  simple 
descriptive  model  for  the  overall  kinematic  behavior  that  captures  these  features  over  most  of  the  probed  volume, 
and  can  be  used  to  search  for  substructure  in  kinematic  and  metallicity  space.  We  use  this  model  to  predict  further 
improvements  in  kinematic  mapping  of  the  Galaxy  expected  from  Gaia  and  the  Large  Synoptic  Survey  Telescope. 

Key  words:  Galaxy:  disk  -  Galaxy:  halo  -  Galaxy:  kinematics  and  dynamics  -  Galaxy:  stellar  content  -  Galaxy: 
structure  -  methods:  data  analysis  -  stars:  statistics 
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1.  INTRODUCTION 


Group  and  mergers  with  neighboring  galaxies.  From  our  van¬ 
tage  point  inside  the  disk  of  the  Milky  Way,  we  have  a  unique 
opportunity  to  study  an  spiral  galaxy  in  great  detail.  By 
measuring  and  analyzing  the  properties  of  large  numbers  of 


The  Milky  Way  is  a  complex  and  dynamic  structure  that  is 
constantly  being  shaped  by  the  infall  of  matter  from  the  Local 
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individual  stars,  we  can  map  the  Milky  Way  in  a  nine¬ 
dimensional  space  spanned  by  the  three  spatial  coordinates, 
three  velocity  components,  and  three  stellar  parameters — lumi¬ 
nosity,  effective  temperature,  and  metallicity. 

In  this  paper,  the  third  in  a  series  of  related  studies,  we  use 
data  obtained  by  the  Sloan  Digital  Sky  Survey  (SDSS;  York 
et  al.  2000)  to  study  in  detail  the  distribution  of  tens  of  millions 
of  stars  in  this  multi-dimensional  space.  In  Juric  et  al.  (2008, 
hereafter  J08),  we  examined  the  spatial  distribution  of  stars  in  the 
Galaxy,  and  in  Ivezic  et  al.  (2008a,  hereafter  108)  we  extended 
our  analysis  to  include  the  metallicity  distribution.  In  this  paper, 
working  with  a  kinematic  data  set  unprecedented  in  size,  we 
investigate  the  distribution  of  stellar  velocities.  Our  data  include 
measurements  from  the  SDSS  astrometric,  photometric,  and 
spectroscopic  surveys:  the  SDSS  Data  Release  7  (Abazajian 
et  al.  2009)  radial-velocity  sample  includes  ~  170,000  main- 
sequence  stars,  while  the  proper-motion  sample  includes  18.8 
million  stars,  with  about  6.8  million  F/G  stars  for  which 
photometric-metallicity  estimates  are  also  available.  These  stars 
sample  a  distance  range  from  ~100  pc  to  ~10  kpc,  probing 
much  farther  from  Earth  than  the  HIPPARCOS  sample,  which 
covers  only  the  nearest  ~100  pc  (e.g.,  Dehnen  &  Binney  1998; 
Nordstrom  et  al.  2004).  With  the  SDSS  data  set,  we  are  offered 
for  the  first  time  an  opportunity  to  examine  in  situ  the  thin/thick 
disk  and  disk/halo  boundaries  over  a  large  solid  angle,  using 
millions  of  stars. 

In  all  three  of  the  papers  in  this  series,  we  have  employed  a 
set  of  photometric-parallax  relations,  enabled  by  accurate  SDSS 
multi-color  measurements,  to  estimate  the  distances  to  main- 
sequence  stars.  With  these  distances,  accurate  to  ~10%-15%, 
the  stellar  distribution  in  the  multi-dimensional  phase  space  can 
be  mapped  and  analyzed  without  any  additional  assumptions. 
The  primary  aim  of  this  paper  is  thus  to  develop  quantitative 
understanding  of  the  large-scale  kinematic  behavior  of  the 
disk  and  halo  stars.  From  the  point  of  view  of  an  observer, 
the  goal  is  to  measure  and  describe  the  radial-velocity  and 
proper-motion  distributions  as  functions  of  the  position  in,  for 
example,  the  r  versus  g  —  r  color-magnitude  diagram,  and  as 
functions  of  the  position  of  the  analyzed  sample  on  the  sky. 
From  the  point  of  view  of  a  theorist,  we  seek  to  directly 
quantify  the  behavior  of  the  probability  distribution  function, 
P(vr,  V(p,  vz\R,  0,  Z,  T ,  L,  [Fe/H]),  where  (vR,  v#,  vz)  are  the 
three  velocity  components  in  a  cylindrical  coordinate  system, 
(R,cj),  Z)  describe  the  position  of  a  star  in  the  Galaxy,  and  T, 
L,  and  [Fe/H]  are  its  temperature,  luminosity,  and  metallicity, 
respectively  (“|”  means  “given”). 

This  a  different  approach  than  that  taken  by  the  widely  used 
“Besan^on”  Galaxy  model  (Robin  &  Creze  1986;  Robin  et  al. 
2003,  and  references  therein),  which  attempts  to  generate  model 
stellar  distributions  from  “first  principles”  (such  as  an  adopted 
initial  mass  function)  and  requires  dynamical  self-consistency. 
Instead,  we  simply  seek  to  describe  the  directly  observed 
distributions  of  kinematic  and  chemical  quantities  without 
imposing  any  additional  constraints.  If  these  distributions  can 
be  described  in  terms  of  simple  functions,  then  one  can  try  to 
understand  and  model  these  simple  abstractions,  rather  than  the 
full  voluminous  data  set. 

As  discussed  in  detail  by  J08  and  108,  the  disk  and  halo 
components  have  spatial  and  metallicity  distributions  that  are 
well  fitted  by  simple  analytic  models  within  the  volume  probed 
by  SDSS  (and  outside  regions  with  strong  substructure,  such  as 
the  Sgr  dwarf  tidal  stream  and  the  Monoceros  stream).  In  this 
paper,  we  develop  analogous  models  that  describe  the  velocity 
distributions  of  disk  and  halo  stars. 


Questions  we  ask  include  the  following: 

1.  What  are  the  limitations  of  the  Schwarzschild  ellipsoidal 
approximation  (a  three-dimensional  Gaussian  distribution, 
Schwarzschild  1979)  for  describing  the  velocity  distribu¬ 
tions? 

2.  Given  the  increased  distance  range  compared  to  older 
data  sets,  can  we  detect  spatial  variation  of  the  best-fit 
Schwarzschild  ellipsoid  parameters,  including  its  orienta¬ 
tion? 

3.  Does  the  halo  rotate  on  average? 

4.  Is  the  kinematic  difference  between  disk  and  halo  stars  as 
remarkable  as  the  difference  in  their  metallicity  distribu¬ 
tions? 

5.  Do  large  spatial  substructures,  which  are  also  traced  in 
metallicity  space,  have  distinctive  kinematic  behavior? 

Of  course,  answers  to  a  number  of  these  questions  are  known 
to  some  extent  (for  excellent  reviews,  see  Gilmore  et  al.  1989; 
Majewski  1993;  Helmi  2008,  for  context  and  references,  see 
also  the  first  two  papers  in  this  series).  For  example,  it  has  been 
known  at  least  since  the  seminal  paper  of  Eggen  et  al.  (1962)  that 
high-metallicity  disk  stars  move  on  nearly  circular  orbits,  while 
many  low-metallicity  halo  stars  move  on  eccentric,  randomly 
oriented  orbits.  However,  given  the  order  of  magnitude  increase 
in  the  number  of  stars  compared  to  previous  work,  larger 
distance  limits,  and  accurate  and  diverse  measurements  obtained 
with  the  same  facility,  the  previous  results  (see  108  for  a 
summary  of  kinematic  results)  can  be  significantly  improved 
and  expanded. 

The  main  sections  of  this  paper  include  a  description  of  the 
data  and  methodology  (Section  2),  followed  by  analysis  of  the 
various  stellar  subsamples.  In  Section  3,  we  begin  by  analyzing 
the  proper-motion  sample  and  determining  the  dependence  of 
the  azimuthal  (rotational)  and  radial-velocity  distributions  on 
position  for  halo  and  disk  subsamples  selected  along  l  =  0° 
and  /  =  180°.  The  spectroscopic  sample  is  used  in  Section  4 
to  obtain  constraints  on  the  behavior  of  the  vertical-velocity 
component,  and  to  measure  the  velocity-ellipsoid  tilts.  The 
resulting  model  is  then  compared  to  the  full  proper-motion 
sample  and  radial- velocity  samples  in  Section  5.  Finally,  in 
Section  6,  we  summarize  and  discuss  our  results,  including  a 
comparison  with  prior  results  and  other  work  based  on  SDSS 
data. 

2.  DATA  AND  METHODOFOGY 

The  characteristics  of  the  SDSS  imaging  and  spectroscopic 
data  relevant  to  this  work  (Fukugita  et  al.  1996;  Gunn  et  al. 
1998, 2006;  Hogg  et  al.  2001;  Smith  et  al.  2002;  Stoughton  et  al. 
2002;  Pier  et  al.  2003;  Ivezic  et  al.  2004;  Tucker  et  al.  2006; 
Abazajian  et  al.  2009;  Yanny  et  al.  2009)  are  described  in  detail 
in  the  first  two  papers  in  the  series  (J08,  108).  Here,  we  only 
briefly  summarize  the  photometric-parallax  and  photometric- 
metallicity  methods,  and  then  describe  the  proper-motion  data 
and  their  error  analysis.  The  subsample  definitions  are  described 
at  the  end  of  this  section. 

2.7.  The  Photometric-parallax  Method 

The  majority  of  stars  in  the  SDSS  imaging  catalogs  (~90%) 
are  on  the  main  sequence  (J08  and  references  therein)  and, 
using  the  broadband  colors  measured  by  SDSS,  it  is  possible  to 
estimate  their  absolute  magnitude.  Briefly,  the  r-band  absolute 
magnitude,  Mr,  of  a  star  can  be  estimated  from  its  position 
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on  the  stellar  locus  of  the  Mr  versus  g  —  i  color-magnitude 
diagram.  The  position  of  this  stellar  locus  is  in  turn  sensitive 
to  metallicity,  so  we  must  apply  an  additional  correction  to  the 
absolute  magnitude.  A  maximum-likelihood  implementation  of 
this  method  was  introduced  and  discussed  in  detail  in  J08. 
The  method  was  further  refined  by  108,  who  calibrated  its 
dependence  on  metallicity  using  globular  clusters. 

We  estimate  absolute  magnitudes  using  Equation  (A7)  in  108, 
which  corrects  for  age  effects,  and  Equation  (A2)  in  the  same 
paper  to  account  for  the  impact  of  metallicity.  The  resulting 
distance  range  covered  by  the  photometric-parallax  relation 
depends  upon  color  and  metallicity,  but  spans  ~100  pc  to 
~10  kpc.  Based  on  an  analysis  of  stars  in  globular  clusters, 
108  estimate  that  the  probable  systematic  errors  in  absolute 
magnitudes  determined  using  these  relations  are  about  0.1  mag, 
corresponding  to  5%  systematic  distance  errors  (in  addition 
to  the  10%-15%  random  distance  errors).  In  addition,  Sesar 
et  al.  (2008,  hereafter  SIJ08)  used  a  large  sample  of  candidate 
wide-binary  stars  to  show  that  the  expected  error  distribution  is 
mildly  non-Gaussian,  with  a  root-mean-square  (rms)  scatter  in 
absolute  magnitude  of  ~0.3  mag.  They  also  quantified  biases  in 
the  derived  absolute  magnitudes  due  to  unresolved  binary  stars. 

2.2.  The  Photometric -metallicity  Method 

Stellar  metallicity  can  significantly  affect  the  position  of  a 
star  in  the  color-magnitude  diagram  (there  is  a  shift  of  ~  1  mag 
between  the  median  halo  metallicity  of  [Fe/H]  ~  —1.5  and 
the  median  disk  metallicity  of  [Fe/H]  ~  —0.2).  SDSS  spec¬ 
troscopy  is  only  available  for  a  small  fraction  of  the  stars  in  our 
sample,  so  we  adopt  a  photometric-metallicity  method  based  on 
SDSS  u  —  g  and  g  —  r  colors.  This  relation  was  originally  cali¬ 
brated  by  108  using  SDSS  spectroscopic  metallicities.  However, 
the  calibration  of  SDSS  spectroscopic  metallicity  changed  at 
the  high-metallicity  end  after  SDSS  Data  Release  6  (Adelman- 
McCarthy  et  al.  2008).  Therefore,  we  recalibrate  their  expres¬ 
sions  as  described  in  the  Appendix.  The  new  calibration,  given  in 
Equation  (Al),  is  applicable  to  F/G  stars  with  0.2  <  g  —  r  <  0.6 
and  has  photometric-metallicity  errors  that  approximately  fol¬ 
low  a  Gaussian  distribution  with  a  width  of  0.26  dex.  In  addition, 
the  ~0.1  dex  systematic  uncertainties  in  SDSS  spectroscopic 
metallicity  (Beers  et  al.  2006;  Allende  Prieto  et  al.  2006,  2008; 
Lee  et  al.  2008a)  are  inherited  by  the  photometric-metallicity 
estimator.  We  emphasize  that  photometric-metallicity  estimates 
are  only  robust  in  the  range  —  2  <  [Fe/H]  <  0  (see  the  Appendix 
for  details). 

For  stars  with  g  —  r  >  0.6,  we  assume  a  constant  metallicity 
of  [Fe/H]  =  —0.7,  motivated  by  results  for  the  disk  metallicity 
distribution  presented  in  108  and  the  fact  that  SDSS  data  are  too 
shallow  to  include  a  large  fraction  of  red  halo  stars.  A  slightly 
better  approach  would  be  to  use  the  disk  metallicity  distribution 
from  108  to  solve  for  best-fit  distance  iteratively,  but  the  resulting 
changes  in  the  photometric  distances  are  negligible  compared 
to  other  systematic  errors. 

2.3.  The  SDSS-POSS  Proper-motion  Catalog 

We  take  proper-motion  measurements  from  the  Munn  et  al. 
(2004)  catalog  (distributed  as  a  part  of  the  public  SDSS  data  re¬ 
leases),  which  is  based  on  a  comparison  of  astrometric  measure¬ 
ments  between  SDSS  and  a  collection  of  Schmidt  photographic 
surveys.  Despite  the  sizable  random  and  systematic  astrometric 
errors  in  the  Schmidt  surveys,  the  combination  of  a  long  base¬ 
line  (~50  years  for  the  POSS-I  survey),  and  a  recalibration  of 


the  photographic  data  using  the  positions  of  SDSS  galaxies  (see 
Munn  et  al.  for  details),  results  in  median  random  proper-motion 
errors  (per  component)  of  only  ~3  mas  yr-1  for  r  <  18  and 
~5  mas  yr-1  for  r  <  20.  As  shown  below,  systematic  errors 
are  typically  an  order  of  magnitude  smaller.  At  a  distance  of 
1  kpc,  a  random  error  of  3  mas  yr-1  corresponds  to  a  veloc¬ 
ity  error  of  ~15  km  s-1,  which  is  comparable  to  the  radial- 
velocity  accuracy  delivered  by  the  SDSS  stellar  spectroscopic 
survey  (~5.3  km  s-1  at  g  =  18  and  20  km  s-1  at  g  =  20.3; 
Schlaufman  et  al.  2009).  At  a  distance  of  7  kpc,  a  random  error 
of  3  mas  yr-1  corresponds  to  a  velocity  error  of  100  km  s-1, 
which  still  represents  a  usable  measurement  for  large  samples, 
given  that  systematic  errors  are  much  smaller  (~20  km  s-1  at 
a  distance  of  7  kpc).  The  small  and  well-understood  proper- 
motion  errors,  together  with  the  large  distance  limit  and  sample 
size  (proper-motion  measurements  are  available  for  about  38 
million  stars  with  r  <  20  from  SDSS  Data  Release  7)  make  this 
catalog  an  unprecedented  resource  for  studying  the  kinematics 
of  Milky  Way  stars. 

We  warn  the  reader  that  proper-motion  measurements  made 
publicly  available  prior  to  SDSS  Data  Release  7  are  known  to 
have  significant  systematic  errors.  Here,  we  use  a  revised  set 
of  proper-motion  measurements  (Munn  et  al.  2008),  which  are 
publicly  available  only  since  Data  Release  7.  As  described  in  the 
next  section,  we  can  assess  the  error  properties  of  this  revised 
proper- motion  catalog  using  objects  with  known  zero  proper 
motion — that  is,  distant  quasars. 

2.3.1.  Determination  of  Proper-motion  Errors  Using  Quasars 

All  known  quasars  are  sufficiently  distant  that  their  proper 
motions  are  vanishingly  small  compared  to  the  expected  ran¬ 
dom  and  systematic  errors  in  the  Munn  et  al.  catalog.  The 
large  number  of  spectroscopically  confirmed  SDSS  quasars 
(Schneider  et  al.  2007)  which  were  not  used  in  the  recalibra¬ 
tion  of  POSS  astrometry  can  therefore  be  used  to  derive  robust 
independent  estimates  of  these  errors.  In  SDSS  Data  Release 
7,  there  are  69,  916  quasars  with  14.5  <  r  <  20,  redshifts  in 
the  range  0.5  <  z  <  2.5,  and  available  proper  motions  (see  the 
Appendix  for  the  SQL  query  used  to  select  and  download  the  rel¬ 
evant  data  from  the  SDSS  CAS).  Within  this  sample  of  quasars, 
the  proper  motions  have  a  standard  deviation  of  ~ 3.1  mas  yr-1 
for  each  component  (determined  from  the  interquartile  range), 
with  medians  differing  from  zero  by  less  than  0.2  mas  yr-1. 
The  dependence  of  the  random  error  on  r-band  magnitude  is 
well-described  by 

erM  =  2.7  +  2.0  x  100-4(r“20)  masyr"1  (1) 

fitting  only  to  quasars  in  the  range  15  <  r  <  20.  When  the  mea¬ 
surements  of  each  proper-motion  component  are  normalized  by 
aM,  the  resulting  distribution  is  approximately  Gaussian,  with 
only  ~1.8%  of  the  quasar  sample  deviating  by  more  than  3a 
from  zero  proper  motion.  In  addition  to  their  dependence  on 
magnitude,  the  random  proper-motion  errors  also  depend  on 
position  on  the  sky,  but  the  variation  is  relatively  small  (~20%, 
see  right  panels  in  Figure  1).  Finally,  we  find  that  the  corre¬ 
lation  between  the  errors  in  the  two  components  is  negligible 
compared  to  the  total  random  and  systematic  errors. 

The  median  proper  motion  for  the  full  quasar  sample  is 
~0.2  mas  yr-1 ,  but  the  systematic  errors  can  be  larger  by  a  factor 
of  2-3  in  small  sky  patches,  as  illustrated  in  Figure  1 .  We  find  that 
the  distribution  of  systematic  proper-motion  errors  in  ~  100  deg2 
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SOSS  0R7  quosors  (/x,  medions:  -1.5  to  1.5  mos/yr)  SDSS  DR7  quosors  (rms  for  /*,:  2  to  5  mos/yr) 


7TCT 


SOSS  0R7  quosors  (/x*  medians:  -1.5  to  1.5  mos/yr) 


SDSS  DR7  quosors  (rms  for  /^:  2  to  5  mos/yr) 


Figure  1.  Behavior  of  proper-motion  measurements  for  60,000  spectroscopically  confirmed  SDSS  quasars  with  b  >  0°.  The  color-coded  maps  (see  the  legend  on 
top,  units  are  mas  yr“ 1 )  show  the  distribution  of  the  median  (left)  and  rms  (right)  for  the  longitudinal  (top)  and  latitudinal  (bottom)  proper-motion  components  in  a 
Lambert  projection  of  the  northern  Galactic  cap.  The  median  number  of  quasars  per  pixel  is  ~250.  For  both  components,  the  scatter  across  the  sky  is  0.60  mas  yr-1. 
The  median  proper  motion  for  the  full  quasar  sample  is  0.15  mas  yr-1  in  the  longitudinal  direction,  and  —0.20  mas  yr-1  in  the  latitudinal  direction.  The  thick  line  in 
the  top-left  panel  shows  the  selection  boundary  for  the  “meridional  plane”  sample. 


patches  of  sky  has  a  width  of  ~0.67  mas  yr- 1  in  each  component, 
about  twice  as  large  as  that  expected  from  purely  statistical 
noise  (per  bin,  using  Equation  (1)).  As  the  figure  shows,  a  few 
regions  of  the  sky  have  coherent  systematic  errors  at  a  level 
~1  mas  yr-1  (e.g.,  the  median  /x/  toward  /  ~  270°,  or  /x&  toward 
the  inner  Galaxy).  Therefore,  the  kinematics  measured  using 
proper  motions  in  these  regions  should  be  treated  with  caution. 

The  largest  systematic  errors,  ~1  mas  yr-1  for  /x/,  are  seen 
toward  l  ~  270°  in  the  top-left  panel  in  Figure  1,  which 
corresponds  to  8  <  10°.  In  this  region,  the  systematic  deviation 
of  quasar  proper  motions  from  zero  is  approximately  parallel 
to  lines  of  constant  right  ascension,  suggesting  that  the  data 
may  be  suffering  from  systematic  effects  due  to  atmospheric 
refraction  and  spectral  differences  between  quasars  and  galaxies 
used  in  the  recalibration  of  POSS  astrometry.  This  effect  would 
be  strongest  for  observations  obtained  at  high  air  mass,  as  are 
typical  for  fields  at  low  declination  (the  POSS  data  were  obtained 
at  a  latitude  of  +33°).  We  find  that  the  median  quasar  proper 
motion  in  the  8  direction  is  well-described  by 

ills)  =  —0.72  +  0.019  5  masyr-1  (2) 


for  —5°  <  8  <  30°,  where  8  is  in  degrees.  At  8  >  30°,  we  find 
<At*>  <  0.2 masyr-1. 

The  observed  direction  and  magnitude  of  this  systematic 
offset  (corresponding  to  an  astrometric  displacement  of  up  to 
~30  mas)  are  consistent  with  detailed  studies  of  atmospheric 
dispersion  effects  on  observations  of  quasars  (Kaczmarczik 
et  al.  2009).  Therefore,  it  is  possible  that  the  systematic  errors 
in  stellar  proper  motions  (whose  spectral  energy  distributions 
differ  less  from  galaxy  spectral  energy  distributions  than  is 
the  case  for  quasars)  are  smaller  than  implied  by  Figure  1. 
Nevertheless,  we  will  conservatively  adopt  the  quasar  proper- 
motion  distributions  as  independent  estimates  of  systematic  and 
random  proper-motion  errors  for  stars  analyzed  in  this  work. 
In  particular,  we  adopt  0.6  mas  yr-1  as  an  estimate  for  typical 
systematic  proper-motion  error. 

The  quasar  sample  has  a  much  narrower  color  range  than 
that  seen  in  main-sequence  stars  (96%  of  the  quasar  sample 
satisfies  —0.2  <  g  —  r  <  0.6),  and  provides  a  better  estimate 
of  systematic  proper-motion  errors  for  the  blue  stars  than  for 
the  red  stars.  Within  the  above  well-sampled  color  range,  we 
find  a  median  proper-motion  gradient  with  respect  to  the  g  —  r 
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color  of  <0.1  mas  yr-1  mag-1  (per  component).  When  the  fit 
is  extended  to  g  —  r  <  1.6  (using  a  much  smaller  number  of 
quasars),  the  gradient  is  still  smaller  than  0.5  mas  yr-1  mag-1. 
Hence,  the  proper-motion  systematics  have  a  color  dependence 
that  is  smaller  than,  or  at  most  comparable  to,  their  dependence 
on  sky  position. 

A  systematic  error  in  proper  motion  of  0.6  mas  yr-1  cor¬ 
responds  to  a  systematic  velocity  error  of  3  km  s-1  at  1  kpc, 
and  ~20  km  s-1  at  7  kpc.  In  addition,  the  ~5%  systematic  dis¬ 
tance  errors  discussed  in  Section  2.1  are  responsible  for  a  ~5% 
systematic  velocity  uncertainty.  Hence,  for  a  disk-like  heliocen¬ 
tric  tangential  velocity  of  20  km  s-1,  proper-motion  systematics 
dominate  at  distances  beyond  ~1  kpc.  Similarly,  for  a  halo-like 
heliocentric  tangential  velocity  of  200  km  s-1,  proper-motion 
systematics  will  dominate  at  distances  greater  than  7  kpc.  At 
smaller  distances,  the  dominant  systematic  in  our  tangential- 
velocity  estimates  comes  from  systematic  distance  errors.  For 
most  of  the  Galaxy  volume  analyzed  in  this  work,  the  systematic 
distance  errors  dominate  over  systematic  proper-motion  errors. 

2.4.  Comparison  of  Proper  Motions  with  Independent 
Measurements 

As  further  tests  of  the  proper-motion  errors,  we  have  analyzed 
two  independent  sets  of  measurements.  As  shown  below,  they 
confirm  the  results  based  on  our  analysis  of  the  quasar  sample. 

We  have  compared  the  SDSS-POSS  proper  motions  to 
proper-motion  measurements  by  Majewski  (1992)  for  a  sam¬ 
ple  of  326  stars  observed  toward  the  North  Galactic  Pole.  The 
measurements  in  the  Majewski  sample  have  random  errors  that 
are  three  times  smaller,  and  comparable,  but  most  likely,  dif¬ 
ferent  systematic  errors.  The  median  proper-motion  differences 
between  the  two  data  sets  are  below  1  mas  yr-1,  with  an  rms 
scatter  3-4  mas  yr-1  (per  coordinate).  Hence,  this  comparison 
is  consistent  with  our  error  estimates  discussed  in  the  preceding 
section,  and  with  the  estimates  of  Majewski  (1992). 

We  have  also  compared  the  SDSS-POSS  proper  motions  with 
proper  motions  from  the  SDSS  stripe  82  region.  In  Bramich 
et  al.  (2008),  proper  motions  are  computed  using  only  SDSS 
data,  and  thus  they  are  expected  to  have  different,  and  probably 
smaller,  systematic  errors  than  the  SDSS-POSS  proper  motions 
(random  errors  for  the  stripe  82  proper  motions  are  larger  by 
about  a  factor  of  2).  For  ~500,000  stars  with  both  SDSS-POSS 
and  Bramich  et  al.  proper-motion  measurements,  we  find  the 
median  differences  and  the  rms  scatter  to  agree  with  expectation. 
A  single  worrisome  result  is  that  the  median  difference  between 
the  two  data  sets  is  a  function  of  magnitude:  we  find  a  gradient  of 
0.8  mas  yr-1  between  r  —  15  and  r  =  20.  It  is  more  likely  that 
this  gradient  is  due  to  systematic  errors  in  centroiding  sources  on 
photographic  plates,  rather  than  a  problem  with  SDSS  data.  This 
gradient  corresponds  to  a  systematic  velocity  error  as  a  function 
of  distance,  Av  ~  4  (D/ kpc)  km  s-1.  For  example,  a  halo  star 
at  5  kpc,  with  a  relative  velocity  of  200  km  s-1,  would  have  a 
systematic  velocity  uncertainty  of  10%.  This  systematic  error 
is  comparable  to  other  sources  of  systematic  errors  discussed 
above,  and  has  to  be  taken  into  account  when  interpreting  our 
results  below. 

2.5.  The  Main  Stellar  Samples 

When  using  proper  motions,  random  errors  in  the  inferred 
velocities  have  a  strong  dependence  on  magnitude,  and  therefore 
distance,  while  systematic  errors  are  a  function  of  position  on 
the  sky,  as  discussed  above.  Random  errors  in  radial-velocity 


measurements  also  depend  on  magnitude,  as  fits  to  spectral 
features  become  more  difficult  at  lower  signal-to-noise  ratios. 
As  such,  when  radial-velocity  and  proper-motion  measurements 
are  analyzed  simultaneously,  the  systematic  and  random  errors 
combine  in  a  complex  way — care  is  needed  when  interpreting 
the  results  of  such  an  analysis. 

In  order  to  minimize  these  difficulties,  we  separately  analyze 
the  proper-motion  sample  and  the  much  smaller  sample  of  stars 
with  radial  velocities.  Furthermore,  motivated  by  the  metallicity 
distribution  functions  (MDFs)  presented  in  108,  we  separately 
treat  the  low-metallicity  “halo”  stars  and  the  high- metallicity 
“disk”  stars.  For  these  two  samples,  we  require  g  —  r  <  0.6, 
the  regime  in  which  the  photometric-metallicity  estimator  is 
believed  to  be  accurate.  Finally,  we  discuss  a  sample  of  “red” 
stars  with  g  —  r  >  0.6  (roughly,  g  —  i  >  0.8),  which  are 
dominated  by  nearby  (<2  kpc)  disk  stars. 

These  samples  are  selected  from  SDSS  Data  Release  7  using 
the  following  common  criteria: 

1.  unique  unresolved  sources  that  show  subarcsecond  paral¬ 
lax:  binary  processing  flags  DEBLENDED_AS_MOVING, 
SATURATED,  BLENDED,  BRIGHT,  and  NODEBLEND 
must  be  false,  and  parameter  nCHILD=0, 

2.  the  interstellar  extinction  in  the  r  band,  Ar  <0.3  mag, 

3.  dust-corrected  magnitudes  in  the  range  14.5  <  r  < 
20  mag, 

4.  high  galactic  latitudes:  \b\  >  20°, 

5.  proper  motion  available, 

yielding  20.1  million  stars.  The  dust  corrections,  Ar,  were 
computed  using  the  Schlegel  et  al.  (1998)  dust  maps,  with 
conversion  coefficients  derived  assuming  an  Rv  =  3.1  dust 
model.  The  intersection  of  the  following  color  criteria  then 
selects  stars  from  the  main  stellar  locus: 

1.  blue  stars  (6.9  million): 

(a)  0.2  <  (g  —  r)  <  0.6, 

(b)  0.7  <  (u—g)  <  2.0 and —0.25  <  (g—r)—0.5(u—g)  < 
0.05, 

(c)  —0.2  <  0.35(g  —  r)  —  (r  —  i)  <  0.10; 

2.  red  stars  (11.9  million): 

(a)  0.6  <  (g  —  r)  <  1.6, 

(b)  -0.15  <  — 0.270r  +  0.800/  -  0.534z  +  0.054  <  0.15, 

where  the  last  condition  is  based  on  a  “principal  color”  or¬ 
thogonal  to  the  stellar  locus  in  the  i  —  z  versus  r  —  i  color-color 
diagram,  as  defined  in  Ivezic  et  al.  (2004).  This  condition  allows 
for  a  0.15  mag  offset  from  the  stellar  locus. 

During  the  analysis,  “blue”  stars  are  often  further  divided  by 
photometric  metallicity  (see  below  for  details)  into  candidate 
“halo”  stars  ([Fe/H]  <  —1.1)  and  candidate  “disk”  stars 
([Fe/H]  >  —0.9).  Subsamples  with  intermediate  metallicities 
include  non-negligible  fractions  of  both  halo  and  disk  stars. 
Although  the  reduced  proper-motion  diagram  is  frequently  used 
for  the  separation  and  analysis  of  these  two  populations,  we  find 
it  inadequate  for  our  purposes;  the  vertical  gradient  in  rotational- 
velocity  blurs  the  kinematic  distinction  between  disk  and  halo 
(for  a  discussion,  see  SIJ08),  and  thus  this  method  is  applicable 
only  to  stars  with  significant  proper  motion  (leading  to  severe 
selection  effects).  Although  metallicities  are  not  available  for 
red  stars,  results  from  108  imply  that  they  are  dominated  by  the 
disk  population  (red  stars  can  only  be  seen  out  to  ~2  kpc). 

For  each  of  the  subsamples  defined  above,  we  further  separate 
those  objects  with  SDSS  spectroscopic  data  (see  the  Appendix 
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for  a  sample  SQL  query)  into  independent  subsamples.  In  total, 
these  spectroscopic  subsamples  include  172,000  stars  (out  of 
352,000  stars  with  spectra),  after  an  additional  requirement  to 
select  only  main- sequence  stars;  that  is,  stars  with  log(g)  > 
3. 5. 27  Of  the  stars  with  spectroscopic  data,  111,000  are  blue 
(0.2  <  g  —  r  <  0.6)  and  61,000  are  red  (0.6  <  g  —  r  <  1.6). 
When  separating  low-  and  high-metallicity  stars  with  spectra, 
we  use  the  spectroscopic  metallicity  (see  Allende  Prieto  et  al. 
2006  for  details).  Due  to  increased  difficulties  with  measuring 
[Fe/H]  for  red  stars  (g  —  r  >  0.6)  from  SDSS  spectra,  we 
adopt  [Fe/H]  =  —0.7  for  all  such  stars;  this  value  is  the  median 
spectroscopic  [Fe/H]  for  stars  with  0.6  <  g  —  r  <  1.3  (cr  = 
0.4  dex).  For  over  90%  of  ~30,000  stars  with  g  —  r  >  1.3, 
[Fe/H]  is  not  successfully  measured. 


(Gunn  et  al.  1979)  and  Hipparcos  measurements  of  Cepheid 
proper  motions  (Feast  &  Whitelock  1997),  we  adopt  vy$R  = 
220  km  s-1  for  the  motion  of  the  local  standard  of  rest  and 
Rq  =  8  kpc  (for  an  analysis  of  other  recent  measurements, 
see  Bovy  et  al.  2009).  The  adopted  value  of  RQ  is  motivated 
by  geometrical  measurements  of  the  motions  of  stars  around 
Sgr  A*,  which  yield  RQ  =  7.94  d=  0.42  kpc  (Eisenhauer  et  al. 
2003).  For  the  solar  peculiar  motion,  we  adopt  the  HIPPARCOS 
result,  i;^’pec  =  —  10.0±0.4kms-1,  i;® ,pec  =  —  5.3±0.6kms-1, 
and  vf,VQC  =  7.2 ±0.4  km  s-1  (Dehnen  &  Binney  1998,  also  see 
Hogg  et  al.  2005).  Using  these  values,  along  with  Equation  (6), 
we  obtain  the  Galactocentric  velocity  components: 

Vi  =  vf*  +  vf,  i  =  X,  Y,  Z,  (7) 


2.6.  Coordinate  Systems  and  Transformations 

Following  J08  and  108,  we  use  a  right-handed,  Cartesian 
Galactocentric  coordinate  system  defined  by  the  following  set 
of  coordinate  transformations: 

X  =  Rq  —  D  cos (/)  cos (b) 

Y  =  —D  sin(/)  cos (b)  (3) 

Z  =  D  sin (/?), 


where  Rq  =  8  kpc  is  the  adopted  distance  to  the  Galactic 
center,  D  is  distance  of  the  star  from  the  Sun,  and  (/ ,  b)  are 
the  Galactic  coordinates.  Note  that  the  Z  =  0  plane  passes 
through  the  Sun,  not  the  Galactic  center  (see  J08),  the  X-axis 
is  oriented  toward  /  =  180°,  and  the  F-axis  is  oriented  toward 
/  =  270°  (the  disk  rotates  toward  /  ~  90°).  The  main  reason 
for  adopting  a  Galactocentric  coordinate  system,  rather  than  a 
traditional  heliocentric  system,  is  that  new  data  sets  extend  far 
beyond  the  solar  neighborhood. 

We  also  employ  a  cylindrical  coordinate  system  defined  by 


R  =  Vx2  +  F2, 


0  =  tan 


(4) 


The  tangential  velocity,  v,  is  obtained  from  the  proper  motion, 
/x,  and  the  distance  D  by 


v  =  4.74  - -  - km  s  l.  (5) 

masyr-1  kpc 

Given  the  line-of-sight  radial  velocity,  ura(j,  and  the  two 
components  of  tangential  velocity  aligned  with  the  Galactic 
coordinate  system,  i p  and  Vb,  the  observed  heliocentric  Cartesian 
velocity  components  are  given  by 

i/^bs  =  —  urad  cos(/)  cos(Z?)  +  Vb  cos(/)  sin(Z?)  +  vi  sin(/) 
u°bs  =  —  t>rad  sin(/)  cos (b)  +  i %  sin(/)  sin(Z?)  —  u/  cos(/)  (6) 
u°bs  =  ymd  sin(Z?)  +  Vb  cos (b). 

These  components  are  related  to  the  traditional  UVW  nomen¬ 
clature  by,  vx  =  —U,vY  =  — V ,  and  vz  =  W9  e.g.,  Binney  & 
Merrifield  (1998). 

In  order  to  obtain  the  Galactocentric  cylindrical  velocity 
components,  we  must  first  correct  for  the  solar  motion.  Taking 
into  account  Hi  measurements  of  the  Galactic  rotation  curve 


where  v%  =  —10  km  s  1,  vf  =  —225  km  s  l,  and  v f  = 
7  km  s-1  (note  that  vf  =  —  ulsr  +  Uy’1560).  Below,  we  dis¬ 
cuss  attempts  to  directly  determine  the  solar  peculiar  motion 
(Section  4.2)  and  ulsr  (Section  5.3)  from  our  data. 

Finally,  the  cylindrical  components,  vR  and  v#,  can  be 
computed  using  a  simple  coordinate  rotation, 


X  Y 
Vr  =  Vx - 1 -  Vy  — 

R  R 

Y  X 

Vfk  ~  —VX  —  +  Vy  — . 
cp  X  R  R 


(8) 


Note  that,  in  our  adopted  system,  the  disk  has  a  prograde  rotation 
Vfp  =  —220  km  s-1;  retrograde  rotation  is  indicated  by  v#  >  0. 
Stars  with  vr  >  0  move  away  from  the  Galactic  center,  and 
stars  with  vz  >  0  move  toward  the  North  Galactic  Pole. 


2. 7.  Analysis  Philosophy 

Such  a  massive  data  set,  extending  to  a  large  distance  limit 
and  probing  a  large  fraction  of  the  Galaxy  volume,  can  be  used 
to  map  the  kinematics  of  stars  in  great  detail.  It  can  also  be 
used  to  obtain  best-fit  parameters  of  an  appropriate  kinematic 
model.  However,  it  is  not  obvious  what  model  (functional  form) 
to  chose  without  at  least  some  preliminary  analysis.  Hence,  in 
the  next  two  sections,  we  first  discuss  various  projections  of  the 
multi-dimensional  space  of  the  available  observable  quantities 
and  obtain  a  number  of  constraints  on  the  spatial  variation  of 
stellar  kinematics.  We  then  synthesize  all  of  the  constraints  into 
a  model  described  in  Section  5.  Before  proceeding  with  our 
analysis,  we  provide  a  brief  summary  of  the  first  two  papers  in 
this  series,  whose  results  inform  our  subsequent  analysis. 

2.8.  A  Summary  of  Papers  I  and  II 

Using  photometric  data  for  50  million  stars  from  SDSS  Data 
Release  4  (Adelman-McCarthy  et  al.  2006),  sampled  over  a 
distance  range  from  100  pc  to  15  kpc,  J08  showed  that  the  stellar 
number  density  distribution,  p(R,  Z,  0)  can  be  well-described 
(apart  from  local  overdensities;  the  J08  best  fit  was  refined  using 
residual  minimization  algorithms)  as  a  sum  of  two  cylindrically 
symmetric  components: 

P(R,  Z,  0)  =  Pd(R ,  Z)  +  pH(R,  Z).  (9) 

The  disk  component  can  be  modeled  as  a  sum  of  two 
exponential  disks 


27  Note  that  the  majority  of  stars  with  g  —  r  >  1.2  do  not  have  reliable 
measurements  of  log(g) — we  assume  that  all  stars  with  g  —  r  >  1.2  are 
main-sequence  stars. 


/  |Z+ZQ|  ( R-Rq )  |Z+Zq  |  (R-Rq)\ 

Pd(R,  Z)  =  pd(Rq)  (e  +eDe  H 2  ^  j , 

(10) 
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while  the  halo  component  requires  an  oblate  power-law  model 

/  R 2  \nH/2 

PH(R,Z)  =  PD{RQ)tH  (r2  +  {. ^ —)  .  (11) 

The  best-fit  parameters  are  discussed  in  detail  by  J08.  We 
have  adopted  the  following  values  for  the  parameters  relevant 
to  this  work  (second  column  in  Table  10  from  J08):  ZQ  = 
25  pc,  H\  =  245  pc,  H2  =  743  pc,  cD  =  0.13,  cH  =  0.0051, 
qH  =0.64,  and  nH  =2.77.  The  normalization  pD(RQ)  (which 
is  essentially  the  local  luminosity  function  for  main-sequence 
stars)  is  listed  in  J08  as  a  function  of  color. 

Using  a  photometric-metallicity  estimator  for  F/G  main- 
sequence  stars,  108  obtained  an  unbiased,  three-dimensional 
metallicity  distribution  of  ~2.5  million  F/G  stars  at  heliocentric 
distances  of  up  to  ~  8  kpc.  They  found  that  the  MDFs  of  the  halo 
and  disk  stars  are  clearly  distinct.  The  median  metallicity  of  the 
disk  exhibits  a  vertical  (with  respect  to  the  Galactic  plane,  Z) 
gradient,  and  no  gradient  in  the  radial  direction  (for  Z  >  0.5  kpc 
and  6  <  R  <  10  kpc). 

Similarly  to  the  stellar  number  density  distribution,  p(R,  Z), 
the  overall  behavior  of  the  MDF  /?([Fe/H]|R,  Z)  for  disk  stars 
can  be  well-described  as  a  sum  of  two  components 

p(x  =  [Fe/H]|R,  Z,  0)  =  [1  -  fH(R,  Z)]  pD(x\Z) 

+  fH(R,  Z)pH(x),  (12) 


where  the  halo  star-count  ratio  is  simply, 


MR,  z)  = 


ph(R,  Z) 

Pd(R ,  Z)  +  Ph(R ,  Z) 


(13) 


The  halo  metallicity  distribution,  pH([ Fe/H]),  is  spatially 
invariant  within  the  probed  volume,  and  well-described  by 
a  Gaussian  distribution  centered  on  [Fe/H]  =  —1.46,  with 
an  intrinsic  (corrected  for  measurement  errors)  width  oH  — 
0.30  dex.  For  \Z\  <  10  kpc,  an  upper  limit  on  the  halo  radial 
metallicity  gradient  is  0.005  dex  kpc-1. 

The  disk  metallicity  distribution  varies  with  Z  such  that  its 
shape  remains  fixed,  while  its  median,  pD,  varies  as 


i  ,  |Z| 

[iD{Z)  =  p oo+A^e  ^  , 


(14) 


with  the  best-fit  parameter  values  H ^  =  0.5  kpc,  =  —0.82, 
and  Am  =  0.55.  The  shape  of  the  disk  metallicity  distribution 
can  be  modeled  as 


pD(x  =  [Fe/H] | Z)  =  0.63  G[x |/x  =  a(Z),  a  =  0.2] 

+  0.37  G[x\p  =  a(Z)  +  0.14,  a  =  0.2], 

(15) 

where  the  position  a  and  the  median  pD  are  related  via 
a(Z)  =  Pd(Z)  —  0.067  (unless  measurement  errors  are  very 
large). 

The  main  result  of  this  third  paper  in  the  series  is  the  extension 
of  these  results  for  number  density  and  metallicity  distributions 
to  include  kinematic  quantities. 

3.  ANALYSIS  OF  THE  PROPER-MOTION  SAMPLE 

We  begin  by  analyzing  the  proper-motion  measurements  of 
stars  observed  toward  the  North  Galactic  Pole.  In  this  region, 
the  Galactocentric  azimuthal  velocity,  v#,  and  radial  velocity, 


vr,  can  be  determined  directly  from  the  proper-motion  mea¬ 
surements  (that  is,  without  knowledge  of  the  spectroscopically 
determined  radial  velocity,  i;rad).  In  this  way,  we  can  study  the 
kinematic  behavior  of  stars  as  a  function  of  metallicity  and  dis¬ 
tance  from  the  Galactic  plane,  Z.  We  then  extend  our  analysis 
to  the  entire  meridional  Y  =  0  plane,  and  study  the  variation 
of  stellar  kinematics  with  R  and  Z.  In  the  following  section,  we 
only  consider  the  northern  Galactic  hemisphere,  where  most  of 
the  proper-motion  data  are  available. 


3.1.  Kinematics  Toward  the  North  Galactic  Pole 


We  select  three  stellar  subsamples  in  the  region  b  >  80°, 
including  14,000  blue  disk  stars  at  Z  <  7  kpc,  23,000  blue 
halo  stars  at  Z  <  7  kpc,  and  a  sample  of  105,000  red  stars  at 
Z  <  1  kpc.  In  Figure  2,  we  plot  the  distribution  of  v#  versus 
vR  for  ^6000  blue  disk  and  halo  stars  at  Z  =  4-5  kpc.  In  this 
and  all  subsequent  two-dimensional  projections  of  the  velocity 
distribution  we  plot  smoothed,  color-coded  maps,  where  the 
velocity  distributions  are  estimated  using  the  Bayesian  density 
estimator  of  Ivezic  et  al.  (2005,  see  their  Appendix  for  the 
derivation  and  a  discussion).  At  a  given  position,  the  density  is 
evaluated  as 


P  = 


(16) 


where  dt  is  the  distance  in  the  velocity-velocity  plane,  and  we 
sum  over  the  A  =  10  nearest  neighbors.  The  normalization 
constant,  C,  is  easily  evaluated  by  requiring  that  the  density 
summed  over  all  pixels  is  equal  to  the  total  number  of  data  points 
divided  by  the  total  area.  The  grid  size  is  arbitrary,  but  the  map 
resolution  is  determined  by  the  density  of  points — we  choose 
pixel  size  equal  to  one  half  of  the  mean  velocity  error.  As  shown 
by  Ivezic  et  al.  (2005),  this  method  is  superior  to  simple  Gaussian 
smoothing.  For  comparison,  we  also  plot  linearly  spaced  density 
contours. 

The  six  panels  of  Figure  2  demonstrate  the  variation  of 
kinematics  with  metallicity,  with  the  full  range  of  metallicities 
(—3  <  [Fe/H]  <  0)  plotted  in  the  upper  left  panel  and 
subsamples  with  increasing  metallicity  running  from  left  to 
right,  top  to  bottom.  The  mean  azimuthal  velocity  varies  strongly 
with  metallicity,  from  a  non-rotating  low-metallicity  subsample 
with  large  velocity  dispersion  (top  center  panel)  to  a  rotating 
high-metallicity  sample  with  much  smaller  dispersion  (bottom- 
right  panel).  This  strong  metallicity-kinematic  correlation  is 
qualitatively  the  same  as  discussed  in  the  seminal  paper  by 
Eggen  et  al.  (1962),  except  that  here  it  is  reproduced  in  situ  with 
a  ~100  times  larger,  nearly  complete  sample,  thus  extending  it 
beyond  the  solar  neighborhood.  There  are  some  indications  of 
substructure  in  the  velocity  distribution,  but  much  of  it  remains 
unresolved  due  to  the  large  velocity-measurement  errors. 

The  substructure  becomes  more  apparent  in  Figure  3,  where 
we  plot  the  same  velocity-space  projection  for  60,000  stars 
within  Z  <2.5  kpc.  In  this  figure,  the  panels  show  subsam¬ 
ples  of  increasing  distance  from  the  Galactic  plane,  beginning 
with  Z  =  0. 1-0.2  kpc  in  the  upper  left  panel  (note  the  chang¬ 
ing  axes  between  the  top  and  bottom  rows).  The  substructure 
seen  in  the  closest  bin  probed  by  red  stars  is  very  similar  to 
the  substructure  seen  in  the  local  HIPPARCOS  sample  (Dehnen 
1998).  These  results  were  based  on  a  maximum-likelihood  anal¬ 
ysis  over  the  entire  sky,  while  our  result  arises  from  a  direct 
mapping  of  the  velocity  distribution  of  stars  selected  from 
a  small  region  (~300  deg2).  Using  a  subsample  of  ~  17,000 
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Figure  2.  Change  of  the  d<j>  vs.  vr  velocity  distribution  with  metallicity,  at  an  approximately  constant  R  and  Z.  Velocities  are  determined  from  proper-motion 
measurements.  The  top-left  panel  shows  the  do  vs.  vr  diagram  for  ~6000  blue  (0.2  <  g  —  r  <  0.4)  stars  at  Z  =  4-5  kpc  and  detected  toward  the  North  Galactic  Pole 
(, b  >  80°).  The  distribution  is  shown  using  linearly  spaced  contours,  and  with  color-coded  maps  showing  smoothed  counts  in  pixels  (low  to  high  from  blue  to  red). 
The  other  five  panels  are  analogous,  and  show  subsamples  selected  by  metallicity,  with  the  [Fe/H]  range  listed  above  each  panel  (also  listed  are  the  median  r-band 
magnitude  and  subsample  size).  The  measurement  errors  are  typically  70  km  s_1  (per  star).  Note  the  strong  variation  of  median  do  with  metallicity. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


HIPPARCOS  stars  with  full  three-dimensional  velocity  infor¬ 
mation,  Nordstrom  et  al.  (2004),  Famaey  et  al.  (2005),  and 
Holmberg  et  al.  (2007,  2009)  have  detected  the  same  kinematic 
morphology.  The  similarity  between  these  HIP  PARC  OS -based 
velocity  distributions  and  ours,  including  the  multi-modal  be¬ 
havior  reminiscent  of  moving  groups  (Eggen  1996),  is  quite 
encouraging,  given  the  vastly  different  data  sources.  The  simi¬ 
larity  of  observed  substructure  with  moving  groups  is  even  more 
striking  for  stars  from  a  closer  distance  bin  (Z  =  50-100  pc), 
matched  to  distances  probed  by  the  HIPPARCOS  sample  (see 
Figure  4).  As  suggested  by  De  Simone  et  al.  (2004),  these  mov¬ 
ing  groups  may  arise  from  irregularities  in  the  Galactic  gravita¬ 
tional  potential. 

The  remainder  of  our  analysis  will  focus  on  blue  stars,  which 
sample  a  much  larger  distance  range.  For  a  detailed  study 
of  the  velocity  distribution  of  nearby  red  stars,  including  a 
discussion  of  non-Gaussianity,  vertex  deviations,  and  difficulties 
with  traditional  thin/thick-disk  separation,  we  refer  the  reader 
to  A.  Kowalski  et  al.  (2010,  in  preparation). 

The  dependence  of  the  rotational  velocity  on  height  above  the 
Galactic  plane  is  shown  in  Figure  5.  The  two  subsamples  display 
remarkably  different  kinematic  behavior  (first  seen  locally  by 
Eggen  et  al.  1962)  with  halo  stars  exhibiting  a  small  constant 
rotational  motion  (~— 20  km  s-1),  and  disk  stars  exhibiting  a 
large  rotational- velocity  component  (~— 200  km  s-1  at  Z  ~ 
1  kpc)  that  decreases  with  height  above  the  Galactic  plane. 

We  have  performed  the  same  analysis  using  proper  motions 
based  only  on  POSS  data,  with  SDSS  positions  not  included  in 


the  proper-motion  fit  (not  publicly  available28).  While  random 
proper-motion  errors  become  larger  when  SDSS  data  are  not 
used,  the  median  rotational  velocity  for  halo  stars  decreases  to 
only  5  km  s-1,  suggesting  that  the  apparent  rotational  motion 
in  the  halo  subsample  is  influenced  by  systematic  errors.  These 
tests  also  suggest  that  the  leading  contribution  to  systematic 
proper-motion  errors  could  be  a  difference  between  the  SDSS 
(digital  data)  and  POSS  (digitized  photographic  data)  centroid- 
determination  algorithms.  In  addition,  Smith  et  al.  (2009)  did 
not  detect  halo  rotation  using  a  smaller  sample,  but  with  more 
robust  proper-motion  measurements  based  on  only  SDSS  data, 
while  Allende  Prieto  et  al.  (2006)  found  no  evidence  for  halo 
rotation  using  SDSS  DR3  radial  velocities.  We  conclude  that 
the  net  halo  rotation  in  the  direction  of  the  North  Galactic  Pole 
is  |  vrot  |  <10  km  s-1.  In  addition,  the  measured  halo  velocity 
dispersion  increases  with  Z,  but  when  random  measurement 
errors  are  taken  into  account,  the  data  are  consistent  with  a 
constant  dispersion  of  a £  =  85  zb  5  km  s-1  (derived  using  the 
test  described  in  Section  5). 

The  decrease  of  rotational  velocity  with  Z  for  disk  stars  (often 
referred  to  as  asymmetric  drift,  velocity  lag,  or  velocity  shear; 
see  Section  3.4  of  108  for  more  details  and  references  to  related 
work)  is  in  agreement  with  a  preliminary  analysis  presented  in 
108.  We  find  that  the  observed  behavior  in  the  Z  =  1-4  kpc 


28 


Available  from  J.  Munn  on  request. 
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Figure  3.  Similar  to  Figure  2,  except  that  the  vs.  vr  velocity  distribution  is  plotted  for  a  range  of  Z.  The  top  row  shows  the  v®  vs.  vr  diagrams  for  ~60,000  red 
(g  —  r  >  0.6)  stars  at  Z  =  100-700  pc  and  observed  toward  the  North  Galactic  Pole.  Each  panel  corresponds  to  a  narrow  range  in  Z,  given  above  each  panel.  The 
measurement  errors  vary  from  ~3  km  s-1  in  the  closest  bin  to  ~12  km  s-1  in  the  most  distant  bin.  Note  the  complex  multi-modal  substructure  in  the  top-left  panel. 
The  bottom  three  panels  are  analogous,  and  show  the  v <j>  vs.  vr  diagrams  for  ~7000  blue  (0.2  <  g  —  r  <  0.4)  stars  with  high  metallicity  ([Fe/H]  >  —0.9).  The 
measurement  errors  vary  from  ~20  km  s-1  in  the  closest  bin  to  ~35  km  s-1  in  the  most  distant  bin.  Note  that  the  median  v&  approaches  zero  as  Z  increases. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


E 


vR  (km/s) 

Figure  4.  Similar  to  the  top-left  panel  in  Figure  3,  except  that  stars  are  selected 
from  a  distance  bin  that  corresponds  to  HIPPARCOS  sample  (Z  =  50-100  pc). 
The  positions  of  Eggen’s  moving  groups  (Eggen  1996)  are  marked  by  circles, 
according  to  the  legend  in  the  bottom-right  corner.  The  horizontal  line  at 
=  —225  km  s_ 1  corresponds  to  vanishing  heliocentric  motion  in  the 
rotational  direction. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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range  can  be  described  by 


(V(P)  =  -205  +  19.2 


km  s  1 


(17) 


The  measured  rotational-velocity  dispersion  of  disk  stars  in¬ 
creases  with  Z  faster  than  can  be  attributed  to  measurement 
errors.  Using  a  functional  form  a  =  a  +  b\Z\c ,  we  obtain  an 
intrinsic  velocity  dispersion  fit  of 


a"  =  30  +  3.0 


Z 

kpc 


2.0 


km  s' 


-l 


(18) 


This  function  and  the  best-fit  rotational  velocity  for  halo  stars 
are  shown  as  dotted  lines  in  the  bottom-right  panel  of  Figure  5 
(see  Table  1  for  a  summary  of  all  best-fit  parameters).  108  fit  a 
linear  model  to  v#  versus  Z,  but  the  difference  between  this  result 
and  their  Equation  (15)  never  exceeds  5  km  s-1  for  Z  <  3  kpc. 
The  errors  on  the  power-law  exponents  of  Equations  (17) 
and  (18)  are  ~0.1  and  ~0.2,  respectively. 

However,  a  description  of  the  velocity  distribution  based 
solely  on  the  first  and  second  moments  (Equations  (17) 
and  ( 1 8))  does  not  fully  capture  the  detailed  behavior  of  our  data. 
As  already  discussed  by  108,  the  rotational- velocity  distribution 
for  disk  stars  is  strongly  non-Gaussian  (see  their  Figure  13).  It 
can  be  formally  described  by  a  sum  of  two  Gaussians,  with  a 
fixed  normalization  ratio  and  a  fixed  offset  of  their  mean  values 
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Figure  5.  Dependence  of  the  rotational  velocity,  uo,  on  distance  from  the 
Galactic  plane  for  14,000  high-metallicity  ([Fe/H]  >  —0.9;  top-left  panel)  and 
23,000  low-metallicity  ([Fe/H]  <  —1.1,  top  right)  stars  with  b  >  80°.  In  the 
top  two  panels,  individual  stars  are  plotted  as  small  dots,  and  the  medians  in 
bins  of  Z  are  plotted  as  large  symbols.  The  2cr  envelope  around  the  medians 
is  shown  by  dashed  lines.  The  bottom  two  panels  compare  the  medians  (left) 
and  dispersions  (right)  for  the  two  subsamples  shown  in  the  top  panels  and  the 
dashed  lines  in  the  bottom  two  panels  show  predictions  of  a  kinematic  model 
described  in  the  text.  The  dotted  lines  in  the  bottom-right  panel  show  model 
dispersions  without  a  correction  for  measurement  errors  (see  Table  1). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


for  |  Z  |  <  5  kpc, 

Pd(x  =  ityl Z)  =  0.75  G[x\vn(Z),  or] 

+  0.25  G[x\vn(Z )  -  34km  s“\  a2],  (19) 


where 


vn(Z)  =  - 194+19.2 


Z 

kpc 


1.25 


km  s 


(20) 


The  intrinsic  velocity  dispersions,  o\  and  a2,  are  modeled  as 
a  +  b\Z\c,  with  best-fit  parameters  listed  in  Table  1  (see  a £  and 
0%).  Closer  to  the  plane,  in  the  0.1  <  Z  <  2  kpc  range  probed  by 
red  stars,  the  median  rotational  velocity  and  velocity  dispersion 
are  consistent  with  the  extrapolation  of  fits  derived  here  using 
much  more  luminous  blue  stars. 

Figure  6  shows  the  v $  distribution  for  four  bins  in  Z  (anal¬ 
ogous  to  Figure  13  from  108),  overlaying  two-component 
Gaussian  fits  with  the  measurement  errors  and  vn(Z)  as  free 
parameters.  The  mean  velocity  and  velocity  dispersion  exhibit 
~10  km  s-1  variations  relative  to  their  expected  values;  while 
such  deviations  could  be  evidence  of  kinematic  substructure, 
they  are  also  consistent  with  the  plausible  systematic  errors. 
We  conclude  that  Equations  (19)  and  (20)  provide  a  good  de¬ 
scription  of  the  disk  kinematics  for  stars  observed  toward  the 
North  Galactic  Pole,  within  the  limitations  set  by  the  random 
and  systematic  errors  in  our  data  set. 

The  fits  to  the  observed  velocity  distributions  for  halo  and 
disk  stars  are  shown  in  Figure  6  and  demonstrate  that  the  vertical 
gradients  in  median  rotational  velocity  and  velocity  dispersion 


Table  1 

Best-fit  Parameters3  for  the  Disk  Velocity  Distribution13 


Quantity 

a 

b 

c 

V(f>1 

-194 

19.2 

1.25 

cr! 

t 

12 

1.8 

2.0 

34 

1.2 

2.0 

4* 

30 

3.0 

2.0 

ctr 

40 

5.0 

1.5 

25 

4.0 

1.5 

Notes.  The  uncertainties  are  typically  ~10  km  s_1  for  a, 
~20%  for  b  and  0. 1-0.2  for  c. 

a  All  listed  quantities  are  modeled  as  a  +  b\Z\c ,  with  Z 
in  kpc,  and  velocities  in  km  s_1 . 

b  The  D0  distribution  is  non-Gaussian,  and  can  be 
formally  described  by  a  sum  of  two  Gaussians  with  a  fixed 
normalization  ratio  A:  1,  with  /*  =  3.0.  The  mean  value 
for  the  second  Gaussian  has  a  fixed  offset  from  the  first 
Gaussian,  v^2  =  (v#1  —  Aity),  with  Av^  =  34  km  s-1. 
Extrapolation  beyond  Z  >  5  kpc  is  not  reliable.  The 
velocity  dispersion  for  the  second  Gaussian  is  given  by 
cr^.  If  this  non-Gaussianity  is  ignored,  the  dispersion 
is  given  by  crP . 


for  disk  stars  seen  in  Figure  5  are  not  due  to  contamination  by 
halo  stars.  To  quantitatively  assess  the  impact  of  “population 
mixing”  as  a  result  of  the  adopted  metallicity-based  classifica¬ 
tion  ([Fe/H]  <  —1.1  for  “halo”  stars  and  [Fe/H]  >  —0.9  for 
“disk”  stars)  on  our  measurements  of  these  gradients,  we  have 
performed  a  series  of  Monte  Carlo  simulations.  Assuming  that 
the  fits  shown  in  Figure  6  accurately  depict  the  intrinsic  veloc¬ 
ity  distributions,  and  adopting  analogous  fits  for  their  metal- 
licity  distributions  (108,  their  Figure  7),  we  have  estimated  the 
expected  bias  in  median  velocity  and  velocity  dispersion  for 
each  population  as  a  function  of  their  relative  normalization.  As 
shown  in  Figure  6  from  108,  the  fraction  of  halo  stars  increases 
with  distance  from  the  plane,  from  about  0.1  at  Z  =  1  kpc  to 
about  0.9  at  Z  =  5  kpc.  We  find  that  the  median  rotation  veloc¬ 
ity  and  velocity  dispersion  biases  are  <10  km  s-1  for  disk  stars 
at  Z  <  3  kpc,  as  well  as  for  halo  stars  at  Z  >  3.5  kpc.  Further¬ 
more,  the  biases  are  <20  km  s-1  for  Z  <  4  kpc  for  disk  stars 
and  to  Z  >  3  kpc  for  halo  stars.  At  Z  =  3  kpc,  the  contamina¬ 
tion  of  both  disk  and  halo  subsamples  by  the  other  population  is 
typically  ~10%-15%.  This  small  contamination  and  the  use  of 
the  median  (as  opposed  to  the  mean)  and  dispersion  computed 
from  the  interquartile  range  results  are  reasonably  small  biases. 
With  the  adopted  metallicity  cutoffs,  the  sample  contamination 
reaches  50%  at  Z  =  1-1.5  kpc  for  halo  stars  and  at  Z  =  4.5  kpc 
for  disk  stars. 

The  dependence  of  the  Galactocentric  radial  velocity  on  Z  is 
shown  for  halo  and  disk  subsamples  in  Figure  7.  The  median 
values  (bottom-left  panel)  are  consistent  with  zero,  within  the 
plausible  systematic  errors  (10-20  km  s_1),  at  all  Z.  The  intrinsic 
dispersion  for  halo  stars  is  consistent  with  a  constant  value  of 
=  135  zb  5  km  s-1.  For  disk  stars,  the  best-fit  functional 
form  a  =  a  +  b\Z\c  is 


40  +  5 


Z 

kpc 


1.5 


km  s' 


-l 


(21) 


The  cr^/cr^  ratio  has  a  constant  value  of  ~  1.3 5  for  Z  <  1.5  kpc, 
and  decreases  steadily  at  larger  Z  to  ~1  at  Z  ~  4  kpc. 
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Figure  6.  Symbols  with  error  bars  are  the  measured  rotational-velocity  distribution,  v&,  for  stars  with  0.2  <  g  —  r  <  0.4,  b  >  80°,  and  Z  =  0.8-1. 2  kpc  (top  left, 
~1500  stars),  1.5-2.0  kpc  (top  right,  ~4100  stars),  3.0-4. 0  kpc  (bottom  left,  ~6400  stars)  and  5.0-7. 0  kpc  (bottom  right,  ~12,500  stars).  The  red  and  green  curves 
show  the  contribution  of  a  two-component  disk  model  (see  Equations  (19)  and  (20)),  the  blue  curves  show  the  Gaussian  halo  contribution,  and  the  magenta  curves  are 
their  sum.  Note  the  difference  in  the  scale  of  the  y-axis  between  the  top  two  panels. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  7.  Analogous  to  Figure  5,  but  for  the  radial- velocity  component,  vr. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


3.2.  Kinematics  in  the  Meridional  Y  ~  0  Plane 

The  analysis  of  the  rotational-velocity  component  can  be 
extended  to  the  meridional  plane  defined  by  Y  =  0,  for  which 
the  longitudinal  proper  motion  depends  only  on  the  rotational- 
velocity  component  and  the  latitudinal  proper  motion,  Vb,  is  a 
linear  combination  of  radial  and  vertical  components, 

Vb  =  sin(b)vR  +  cos  (b)vz-  (22) 


Figure  8  plots  v#  and  Vb  as  functions  of  R  and  Z  for  halo  and 
disk  stars  within  10°  of  the  meridional  plane.  The  median  Vb  is 
close  to  zero  throughout  most  of  the  plotted  region,  as  would 
be  expected  if  the  median  vr  and  vz  are  zero  (the  behavior 
of  vz  is  discussed  in  the  next  section).  One  exception  is  a 
narrow  feature  with  Vb  ~  —100  km  s-1  for  R  <  4  kpc. 
While  a  cold  stellar  stream  could  produce  such  a  signature,  its 
narrow  geometry  points  directly  at  the  observer.  This  behavior 
is  also  consistent  with  a  localized  systematic  proper-motion 
error.  Indeed,  the  bottom-left  panel  in  Figure  1  shows  that  the 
systematic  latitudinal  proper- motion  error  at  /  ~  0°,  b  ~  45° 
is  about  1  mas  yr-1,  corresponding  to  a  velocity  error  of 
~100  km  s-1  at  a  distance  of  7  kpc. 

As  seen  in  the  upper  left  panel  of  Figure  8,  the  median  v $  for 
halo  stars  is  close  to  zero  for  R  <  12  kpc.  In  the  region  with 
R  >  12  kpc  and  Z  <  6  kpc,  the  median  indicates  a  surprising 
prograde  rotation  in  excess  of  100  km  s-1 .  This  behavior  is  also 
seen  in  disk  stars,  and  is  likely  due  to  the  Monoceros  stream, 
which  has  a  metallicity  intermediate  between  disk  and  halo  stars 
and  rotates  faster  than  disk  stars  (see  Sections  3.5.1  and  3.5.2  in 
108).  There  is  also  an  indication  of  localized  retrograde  rotation 
for  halo  stars  with  Z  ~  9  kpc  and  R  ~  15  kpc  (corresponding 
to  l  ~  180°,  b  ~  50°,  and  a  distance  of  ~11  kpc).  Stars  with 
Z  =  8-10  kpc  and  R  =  15-17  kpc  have  median  larger  by 
40  km  s-1  (a  ~lcr  effect)  and  median  [Fe/H]  larger  by  0.1  dex 
(~5a  effect)  than  stars  with  Z  =  8-10  kpc  and  R  =  7-13  kpc. 
A  systematic  error  in  /z/  of  ~0.8  mas  yr-1  is  required  to  explain 
this  kinematic  feature  as  a  data  problem  (although  this  would 
not  explain  the  metallicity  offset).  However,  the  top-right  panel 
in  Figure  1  shows  that  the  systematic  /z/  errors  in  this  sky  region 
are  below  0.5  mas  yr-1,  so  this  feature  may  well  be  real.  We 
note  that  in  roughly  the  same  sky  region  and  at  roughly  the  same 
distance,  Grillmair  &  Dionatos  (2006)  have  detected  a  narrow 
stellar  stream. 
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Figure  8.  Dependence  of  velocity,  measured  using  proper  motions,  on  cylindrical  Galactocentric  coordinates  for  172,000  metal-poor  halo  stars  ([Fe/H]  <  —1.1; 
top  panels)  and  205,000  metal-rich  disk-like  stars  ([Fe/H]  >  —0.9;  bottom  panels).  Stars  are  selected  from  three  regions:  b  >  80°  (the  North  Galactic 
pole),  170°  <  /  <  190°  (anticenter),  and  350°  <  /  <  10°  (Galactic  center).  The  left  column  plots  rotational  velocity,  i>o,  while  the  right  column  plots 
vb  =  sin (b)vR  +  cos(b)vz .  To  aid  visualization  of  these  boundaries,  see  the  thick  line  in  the  top-left  panel  in  Figure  1.  The  median  values  of  velocity  in  each 
bin  are  color-coded  according  to  the  legend  shown  in  each  panel.  The  measurements  are  reliable  for  distances  up  to  about  7  kpc,  but  regions  beyond  this  limit  are 
shown  for  halo  stars  for  completeness.  The  fraction  of  disk  stars  is  negligible  at  such  distances;  their  velocity  distribution  is  shown  for  Z  <  6  kpc.  The  region  with 
negative  velocity  on  the  right  side  of  top-left  panel  is  due  to  contamination  of  the  halo  sample  by  stars  from  the  Monoceros  stream.  The  thin  region  with  negative 
velocity  on  the  left  side  of  top-right  panel  is  a  data  artifact  (see  the  text). 


In  order  to  visualize  the  extent  of  “contamination”  by  the 
Monoceros  stream,  we  replace  the  rotational  velocity  for  each 
disk  star  by  a  simulated  value  drawn  from  the  distribution 
described  by  Equation  (19).  We  subtract  this  model  from  the 
data,  and  the  residuals  are  shown  in  the  right-hand  panel  of 
Figure  9.  The  position  of  the  largest  deviation  is  in  excellent 
agreement  with  the  position  of  Monoceros  stream  quantified  in 
108  ( R  =  15-16  kpc  and  |Z|  ~  3-5  kpc).  Further  evidence  for 
the  presence  of  the  Monoceros  stream  is  shown  in  Figure  10, 
in  which  we  analyze  v $  versus  [Fe/H],  as  a  function  of  R ,  for 
blue  stars  at  Z  =  4-6)  kpc.  As  is  evident  in  the  bottom-right 
panel,  there  is  a  significant  excess  of  stars  at  R  >  17  kpc  with 
—  1.5  <  [Fe/H]  <  —0.5  that  rotate  in  a  prograde  direction  with 
~200  km  s_1. 

4.  ANAFYSIS  OF  THE  SPECTROSCOPIC  SAMPFE 

Despite  its  smaller  size,  the  SDSS  DR7  spectroscopic  sample 
of  ~  100,000  main-sequence  stars  is  invaluable,  because  it 
enables  a  direct29  study  of  the  three-dimensional  velocity 
distribution.  The  sample  extends  to  a  distance  of  ~10  kpc, 
at  which  it  can  deliver  velocity  errors  as  small  ~10  km  s-1 


29  Statistical  deprojection  methods,  such  as  that  recently  applied  to  a 
subsample  of  M  stars  discussed  by  Fuchs  et  al.  (2009),  can  be  used  to 
indirectly  infer  the  three-dimensional  kinematics  from  proper-motion  data. 


(corresponding  tangential- velocity  errors  are  ~150  km  s-1  at  a 
distance  of  10  kpc).  For  each  object  in  the  SDSS  spectroscopic 
survey,  its  spectral  type,  radial  velocity,  and  radial-velocity  error 
are  determined  by  matching  the  measured  spectrum  to  a  set 
of  stellar  templates,  which  were  calibrated  using  the  EFODIE 
stellar  library  (Prugniel  et  al.  2007).  Random  errors  on  the  radial- 
velocity  measurements  are  a  strong  function  of  spectral  type  and 
signal-to-noise  ratio,  but  are  usually  <5  km  s-1  for  stars  brighter 
than  g  ~  18,  rising  sharply  to  ~15  km  s-1  for  stars  with  r  =  20. 
We  model  the  behavior  of  the  radial-velocity  errors  as 

<7rad  =  3  +  12  X  10°-4  (r_20)  km  s-1.  (23) 

We  begin  our  analysis  with  blue  disk  and  halo  stars,  and  then 
briefly  discuss  the  kinematics  of  nearby  red  M  stars. 

4.1.  Velocity  Distributions 

We  select  1 1 1,000  stars  with  0.2  <  g  —  r  <  0.6  (74,000  have 
b  >  20°)  and,  using  their  spectroscopic  metallicity,  separate 
them  into  47,000  candidate  halo  stars  with  [Fe/H]  <  —1.1, 
and  ^53,000  disk  stars  with  [Fe/H]  >  —0.9.  Assuming  the 
spectroscopic  metallicities  accurately  separate  disk  from  halo 
stars,  the  use  of  photometric  metallicity  for  the  same  selection 
would  result  in  a  contamination  rate  of  6%  for  the  halo 
subsample,  and  12%  for  the  disk  subsample. 
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Figure  9.  Left  panel  is  analogous  to  the  bottom-left  panel  in  Figure  8,  but  for  the  model  described  in  the  text.  The  right  panel  shows  the  median  difference  between  the 
data  and  model.  Large  discrepancies  at  R  >  12  kpc  are  due  to  the  Monoceros  stream  (at  R  =  18  kpc  and  Z  =  4  kpc;  disk  stars  rotate  with  a  median  v 0  ~  — 100  km  s_1 , 
while  for  the  Monoceros  stream,  v $  ~  —200  km  s-1). 
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Figure  10.  Distribution  of  stars  with  0.2  <  g  —  r  <  0.4  and  Z  =  4-6  kpc  in  the  rotational  velocity  vs.  metallicity  plane,  for  four  ranges  of  Galactocentric  cylindrical 
radius,  R  (top  left:  3-4  kpc;  top  right:  7-9  kpc;  bottom  left:  12-13  kpc;  bottom  right:  17-19  kpc).  In  each  panel,  the  color-coded  map  shows  the  logarithm  of  counts 
in  each  pixel,  scaled  by  the  total  number  of  stars.  The  horizontal  lines  at  v&  =  0  km  s-1  and  v&  =  —220  km  s-1  are  added  to  guide  the  eye.  High-metallicity 
([Fe/H]  ~  —  1)  stars  with  fast  rotation  (i><j>  ~  —220  km  s-1)  visible  in  the  bottom-right  panel  belong  to  the  Monoceros  stream,  and  are  responsible  for  the  features 
seen  at  R  >  15  kpc  in  the  two  left  panels  in  Figure  8. 
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Figure  11.  Similar  to  Figure  5,  but  for  the  vertical-velocity  component, 
vz,  and  using  a  sample  of  stars  with  SDSS  radial-velocity  measurements, 
0.2  <  g  —  r  <  0.4  and  b  >  20°  (12,000  stars  in  the  high-metallicity 
subsample,  and  38,000  stars  in  the  low-metallicity  subsample).  An  analogous 
figure  for  extended  samples  of  53,000  disk  stars  and  47,000  halo  stars  with 
0.2  <  g  —  r  <  0.6  has  a  similar  appearance.  The  behavior  of  the  rotational- 
and  radial-velocity  components  in  this  sample  is  consistent  with  that  shown  in 
Figures  5  and  7. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


The  dependence  of  the  median  vertical  velocity,  vz,  and  its 
dispersion  on  height  above  the  disk,  is  shown  in  Figure  1 1  for 
the  halo  and  disk  subsamples.  The  median  values  of  vz  are 
consistent  with  zero  to  better  than  10  km  s_1  at  Z  <  5  kpc, 
where  statistical  fluctuations  are  small. 

As  with  cr^  and  aR ,  the  vertical  velocity  dispersion  can  be 
modeled  using  a  constant  dispersion  for  halo  stars  (cr^  = 
85  km  s-1),  while  for  disk  stars,  the  best-fit  functional  form 
is 


cf%  =  25  +  4 


Z  1,5 
kpc 


km  s  l. 


(24) 


The  other  two  velocity  components  behave  in  a  manner 
consistent  with  Equations  (17),  (18),  and  (21),  just  as  they 
did  in  the  proper-motion  sample.  This  is  encouraging,  because 
the  spectroscopic  sample  is  collected  over  the  entire  northern 
hemisphere,  unlike  the  proper-motion  subsample  studied  in 
Section  3.1,  which  is  limited  to  b  >  80°. 

The  availability  of  all  three  velocity  components  in  the  spec¬ 
troscopic  sample  makes  it  possible  to  study  the  orientation  of 
the  halo  velocity  ellipsoid.  Figure  12  shows  two-dimensional 
projections  of  the  velocity  distribution  for  subsamples  of  candi¬ 
date  halo  stars  with  0.2  <  g  —  r  <  0.4.  The  top  row  corresponds 
to  stars  above  the  Galactic  plane  at  3  <  Z /kpc  <  4,  while  the 
bottom  row  is  for  stars  the  same  distance  below  the  plane.  The 
velocity  ellipsoid  is  clearly  tilted  in  the  top-  and  bottom-left 
panels,  with  a  tilt  angle  consistent  with  tan~l(vz/vR)  =  R/z. 
While  the  tilt-angle  errors  are  too  large  to  obtain  an  improve¬ 
ment  over  existing  measurements  of  RQ,  it  is  remarkable  that  the 


northern  and  southern  subsamples  agree  so  well.30  In  addition, 
when  the  Z  =  3-5  kpc  sample  is  divided  into  three  subsamples 
with  7  <  R/ kpc  <11,  the  tilt  angle  varies  by  the  expected  ~8° 
in  the  correct  direction  (see  Figure  13).  For  all  bins  in  the  R-Z 
plane,  the  best-fit  tilt  angle  is  statistically  consistent  (within  5°) 
with  tan-1  (vz/vr)  =  R/z.  The  other  two  projections  of  the 
velocity  distribution  for  halo  stars  do  not  exhibit  significant  tilts 
to  within  ~3°. 

If  we  transform  the  velocities  to  a  spherical  coordinate  system, 


R  Z 

Vr  =  Vr~r~  +  Vz~r~ 

gc  ^gc 


ve 


vr 


R , 


—  Vz 


gc 


R 

Rgc 


(25) 


where  r  =  Rgc  =  ( R 2  +  Z2)1/2  is  the  spherical  Galactocentric 
radius,  we  find  no  statistically  significant  tilt  in  any  of  the  two- 
dimensional  velocity- space  projections  for  halo  stars  (with  tilt- 
angle  errors  ranging  from  ~l°to~5°). 

As  shown  in  Figure  14,  we  see  no  evidence  for  a  velocity- 
ellipsoid  tilt  in  vz  versus  vR  for  the  disk  stars.  The  plotted 
subsamples  are  again  selected  to  have  colors  0.2  <  g  —  r  <  0.4, 
but  are  selected  closer  to  the  Galactic  plane,  \Z\  =  1.5-2. 5  kpc, 
in  order  to  improve  statistics  and  reduce  contamination  from 
halo  stars.  The  velocity-ellipsoid  tilt  is  consistent  with  zero 
within  ~lcr,  and  alignment  of  the  velocity  ellipsoid  with  the 
spherical  coordinate  system  of  Equation  (25)  is  ruled  out  at  a 
~2 o  or  greater  confidence  level  for  each  of  five  analyzed  R-Z 
bins  ( R  =  6-1 1  kpc,  with  AR  =  1  kpc).  We  conclude  that  there 
is  no  evidence  for  a  velocity-ellipsoid  tilt  in  the  disk  subsample, 
but  caution  that,  due  to  the  small  Z  range,  the  data  cannot 
easily  distinguish  between  cylindrical  and  spherical  alignment. 
A  model  prediction  for  velocity-ellipsoid  tilt  is  discussed  below. 

The  vertex  deviation  is  analogous  to  the  velocity-ellipsoid  tilt 
discussed  above,  but  is  defined  in  the  v $  versus  vR  plane  instead 
of  the  vz  versus  vR  plane.  The  same  plots  for  red  (g  —  r  >  0.6, 
median  1.2)  disk  stars  are  shown  in  the  center  top  and  bottom 
panels  in  Figure  15.  These  stars  can  be  traced  closer  to  the  plane, 

|  Z |  =  0.6-0. 8  kpc;  in  both  hemispheres,  the  data  are  consistent 
with  a  vertex  deviation  of  ~20°,  with  an  uncertainty  of  ~10°. 
This  result  is  consistent  with  the  vertex  deviation  obtained  by 
Fuchs  et  al.  (2009). 

Another  interpretation  for  the  v $  versus  vR  distribution  of 
disk  stars  invokes  a  two-component  velocity  distribution,  which 
can  result  in  a  similar  deviation  even  if  each  component 
is  perfectly  symmetric  in  the  cylindrical  coordinate  system. 
A.  Kowalski  et  al.  (2010,  in  preparation)  find  that  the  v $ 
versus  vR  distribution  for  red  stars  toward  the  North  Galactic 
Pole,  with  0.1  <  Z/ kpc  <  1.5,  can  be  fitted  by  a  sum  of 
two  Gaussian  distributions  that  are  offset  from  each  other  by 
~10  km  s-1  in  each  direction.  This  offset  results  in  a  non¬ 
zero  vertex  deviation  if  the  sample  is  not  large  enough,  or  if 
measurements  are  not  sufficiently  accurate,  to  resolve  the  two 
Gaussian  components.  This  double- Gaussian  structure  would  be 
at  odds  with  the  classical  description  based  on  the  Schwarzschild 
approximation — we  refer  the  interested  reader  to  the  Kowalski 
et  al.  study  for  more  details.  Unfortunately,  the  spectroscopic 
samples  are  not  large  enough  to  distinguish  a  two-component 
model  from  the  standard  interpretation. 


30  A  plausible,  if  somewhat  optimistic,  tilt-angle  uncertainty  of  1° 
corresponds  to  an  RQ  error  of  0.5  kpc;  extending  the  sample  to  |Z|  =  8  kpc 
could  deliver  errors  of  0.3  kpc  per  bin  of  a  similar  size. 
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Figure  12.  Three  two-dimensional  projections  of  the  velocity  distribution  for  two  subsamples  of  candidate  halo  stars  selected  using  spectroscopic  metallicity 
(—3  <  [Fe/H]  <  —1.1)  and  with  6  <  /?/kpc  <11.  The  top  row  corresponds  to  2700  stars  with  distances,  3  <  Z/kpc  <  4,  and  the  bottom  row  to  1300  stars  with 
— 4  <  Z /kpc  <  —  3.  The  distributions  are  shown  using  linearly  spaced  contours,  and  with  a  color-coded  map  showing  smoothed  counts  in  pixels  (low  to  high  from  blue 
to  red).  The  measurement  errors  are  typically  60  km  s-1.  Note  the  strong  evidence  for  a  velocity-ellipsoid  tilt  in  the  top  and  bottom-left  panels  (see  also  Figure  13). 
The  two  dashed  lines  in  these  panels  show  the  median  direction  toward  the  Galactic  center. 
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Figure  13.  Illustration  of  the  velocity-ellipsoid  tilt-angle  variation.  Analogous  to  Figure  12,  except  that  only  the  vz  vs.  vr  projection  is  shown  for  a  constant  Z,  for 
three  ranges  of  R,  as  marked  on  the  top  of  each  panel. 


4.1.1.  A  Model  Prediction  for  Velocity -Ellipsoid  Tilt 

The  tilt  of  the  velocity-ellipsoid  tilt  in  the  vz~vr  plane  is 
calculated  using  methods  developed  in  Kent  &  de  Zeeuw  (1991), 
who  show  that  the  tilt  of  the  ellipsoid  at  any  point  in  the  Milky 
Way  depends  not  only  on  the  gravitational  potential  of  the 
Galaxy  but  also  on  the  isolating  integrals  of  a  particular  orbit  and 
thus  (weakly)  on  the  distribution  function  of  stellar  velocities. 
We  use  the  model  of  the  Milky  Way  gravitational  field  from 
Carlberg  &  Innanen  (1987),  scaled  to  a  solar  radius  of  8  kpc 


and  a  circular  velocity  of  220  km  s_1.  At  a  particular  point 
in  the  Galaxy,  we  integrate  eight  orbits,  each  one  launched 
from  that  point,  with  velocities  in  each  of  the  R ,  Z,  and  0 
coordinates  that  are  zb  la  from  the  systemic  velocity.  We  use  the 
observed  velocity  distributions  to  obtain  the  systemic  velocity 
and  the  1  a  values.  Furthermore,  using  the  “least-squares  fitting” 
method,  we  determine  the  parameters  of  the  prolate  spheroidal 
coordinate  system  that  best  matches  each  orbit — from  this,  we 
can  compute  its  contribution  to  the  velocity  ellipsoid.  For  some 
orbits  (especially  those  of  halo  stars),  the  orbits  are  not  local  to 
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Figure  14.  Analogous  to  Figure  12,  except  that  the  velocity  distribution  is  shown  for  two  subsamples  of  candidate  disk  stars  selected  using  spectroscopic  metallicity 
([Fe/H]  >  —0.9).  The  top  row  corresponds  to  2200  stars  with  distances  from  the  Galactic  plane  1.5  <  Z/kpc  <  2.5,  and  the  bottom  row  to  1500  stars  with 
—2.5  <  Z/kpc  <  —1.5.  The  measurement  errors  are  typically  35  km  s_1 .  Note  the  absence  of  velocity-ellipsoid  tilt  in  the  top  and  bottom-left  panels. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  15.  Analogous  to  Figure  14,  except  that  the  velocity  distribution  is  shown  for  two  subsamples  of  red  stars  (g  —  r  >  0.6):  the  top  row  corresponds  to  3000  stars 
with  distances  from  the  Galactic  plane  0.6  <  Z/kpc  <  0.8,  and  the  bottom  row  to  4600  stars  with  —0.8  <  Z/kpc  <  —0.6.  The  measurement  errors  are  typically 
~15  km  s_1 .  There  is  no  strong  evidence  in  these  panels  for  a  tilt  in  the  velocity  ellipsoid. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Table  2 

A  Model  Prediction  for  Velocity-Ellipsoid  Tilta 


R  (kpc) 

Z  (kpc) 

Orz  (deg) 

arctan (Z/R)  (deg) 

8.0 

0.7 

3 

5.0 

8.0 

2.0 

11 

14.0 

8.0 

3.5 

26 

23.6 

7.0 

4.0 

32 

29.7 

9.0 

4.0 

25 

24.0 

11.0 

4.0 

21 

20.0 

Notes.  a  The  first  two  columns  determine  the  position  in 
the  Galaxy.  The  third  column  lists  the  predicted  orien¬ 
tation  of  the  velocity  ellipsoid  computed  as  described  in 
Section  4.1.1.  The  fourth  column  lists  the  orientation  of 
a  velocity  ellipsoid  pointed  toward  the  Galactic  center. 

the  solar  neighborhood,  and  the  local  fitting  method  cannot  fit 
the  entire  orbit.  For  these  cases,  we  confine  the  orbits  to  radii 
greater  than  4  kpc.  Since  the  only  purpose  of  the  fit  is  to  model 
the  kinematics  in  a  small  volume  of  space,  the  fact  that  the 
fit  is  no  longer  global  is  of  no  consequence.  The  scatter  in  tilt 
angle  among  the  individual  orbits  at  each  position  depends  only 
weakly  on  position  (~10%). 

Table  2  gives  the  predicted  tilt  angles  for  each  position  plotted 
in  Figures  12-15.  For  the  halo  stars  plotted  in  Figures  12 
and  13,  where  the  ellipsoid  shows  the  most  distinctive  tilt,  the 
predictions  match  the  observed  tilts  quite  well,  but  are  slightly 
steeper  at  Z  >  3  kpc  than  the  Z/R  relation  predicted  for 
spherical  symmetry.  For  the  disk  stars  in  Figures  14  and  15,  the 
predicted  tilts  lie  between  the  cases  of  cylindrical  and  spherical 
symmetry.  The  observed  ellipsoid  is  sufficiently  round,  however, 
that  no  definitive  comparison  can  be  made. 

4.2.  Direct  Determination  of  the  Solar  Peculiar  Motion 

If  there  is  no  net  streaming  motion  in  the  Z  direction  in  the 
solar  neighborhood,  the  median  heliocentric,  i^bs,  for  nearby 
stars  should  be  equal  to  v®  (7  km  s_1,  based  on  an  analysis 
of  HIPPARCOS  results  by  Dehnen  &  Binney  1998).  We  do 
not  expect  a  large  velocity  gradient  within  ~1  kpc  from  the 
Sun,  so  we  select  from  the  spectroscopic  sample  ~  13,000  M 
dwarfs  with  2. 3  <  g  —  i  <2. 8.  In  the  northern  hemisphere,  we 
have  5700  stars  with  a  median  heliocentric  velocity,  (t5?s)  = 
—  1.8  km  s-1,  while  for  stars  in  the  southern  hemisphere  we 
obtain  (u^bs)  =  —11.0  km  s-1.  This  difference  is  likely  due 
to  a  systematic  radial- velocity  error.  If  we  simultaneously  vary 
an  assumed  radial- velocity  error,  Arad,  and  the  solar  peculiar 
motion,  i>®,  while  requiring  that  the  median  i^bs  should  be  the 
same  for  both  hemispheres,  we  obtain  Arad  =  5.0  zb  0.4  km  s-1 
and  u®  =  6.5  zb  0.4  km  s-1.  This  value  for  i>®  is  in  excellent 
agreement  with  the  HIPPARCOS  value  of  7.2  ±  0.4  km  s-1 
(Dehnen  &  Binney  1998).  This  systematic  offset  in  SDSS  radial 
velocities  is  probably  due  to  the  small  number  of  ELODIE 
templates  for  red  stars.  It  is  likely  that  the  adoption  of  improved 
templates  from  Bochanski  et  al.  (2007a)  will  yield  smaller 
systematic  errors.  A  similar  analysis  for  blue  stars  does  not  yield 
a  robust  detection  of  the  velocity  offset.  A  detailed  comparison 
of  SDSS  radial  velocities  with  radial- velocity  standards  from 
the  literature  arrived  at  the  same  null  result  for  blue  stars. 

As  with  the  vertical  component  of  the  solar  peculiar  motion, 
if  the  adopted  value  of  v%  =  —10  km  s-1  were  incorrect,  the 
median  vr  would  deviate  from  zero.  The  rms  scatter  of  the 
median  vR  for  subsamples  of  nearby  M  stars  selected  by  distance 
and  color  is  0.5  km  s-1,  which  is  an  upper  limit  on  the  error  in 


the  adopted  value  of  u®.  This  result,  which  is  based  on  the  full 
three-dimensional  velocity  distribution,  agrees  well  with  results 
from  indirect  statistical  deprojection  methods  using  only  proper 
motions  (Dehnen  1998;  Fuchs  et  al.  2009). 

For  both  blue  and  red  disk  stars,  the  extrapolation  of  the 
median  to  Z  =  0  yields  —205  km  s-1.  Since  we  corrected 
stellar  velocities  for  an  assumed  solar  motion  of  —225  km  s-1, 
this  implies  that  the  Y  component  of  the  solar  velocity  relative  to 
the  bulk  motion  of  stars  in  the  solar  neighborhood  is  20  km  s_1 , 
in  agreement  with  recent  results  obtained  by  Fuchs  et  al.  (2009) 
for  the  same  data  set.  A  similar  value  was  obtained  by  Dehnen  & 
Binney  (1998)  for  their31  subsample  of  red  stars  within  100  pc. 

5.  A  MODEL  FOR  THE  KINEMATICS  OF  DISK  AND 
HALO  STARS 

Informed  by  the  results  from  the  preceding  two  sections,  we 
introduce  a  model  that  aims  to  describe  the  global  behavior  of  the 
observed  stellar  kinematics.  In  our  model,  we  do  not  attempt  to 
account  for  kinematic  substructure  (e.g.,  the  Monoceros  stream), 
or  the  Galactic  bulge  region,  nor  do  we  incorporate  any  complex 
kinematic  behavior  close  to  the  Galactic  plane.  Nevertheless,  we 
attempt  to  capture  the  gross  properties  of  the  data  in  the  volume 
probed  by  SDSS,  including  the  bulk  kinematic  trends  and  the 
kinematic  differences  between  high-metallicity  disk  stars  and 
low-metallicity  halo  stars.  We  describe  the  model  in  Section  5.1, 
then  test  it  in  Section  5.2  using  both  the  proper-motion  and 
radial- velocity  samples. 

5.7.  The  Kinematic  Model 

For  halo  stars,  a  single  velocity  ellipsoid  (expressed  in  a 
Galactocentric  spherical  coordinate  system,  see  Section  4.1)  is  a 
good  description  of  the  gross  halo  kinematics  within  the  10  kpc 
distance  limit  of  our  sample.  Our  model  assumes  that  the  halo 
has  no  net  rotation  (see  below  for  a  test  of  this  assumption),  and 
that  the  principal  axes  of  the  velocity  ellipsoid  are  aligned  with  a 
spherical  coordinate  system.  The  velocity  dispersions  measured 
in  Sections  3  and  4,  a £  =  135  km  s_1,  =  85  km  s_1, 

and  cr^  =  85  km  s_1,  are  expressed  in  a  cylindrical  coordinate 
system.  The  interplay  between  proper-motion,  radial-velocity, 
and  distance  measurement  errors  is  complex,  so  we  use  Monte 
Carlo  simulations  to  translate  them  to  spherical  coordinates; 
the  results  of  this  exercise  are  atH  =  141  km  s_1  and  <Jq  = 
75  km  s_1  ( o £  is  unchanged),  with  uncertainties  of  ~5  km  s_1. 

For  disk  stars  within  ~l-2  kpc  from  the  Sun,  the  velocity- 
measurement  errors  are  sufficiently  small,  and  the  samples  are 
sufficiently  large,  to  resolve  rich  kinematic  substructure  (e.g., 
Figure  3).  This  behavior  is  quantified  in  detail  in  A.  Kowalski 
et  al.  (2010,  in  preparation),  but  here  we  simply  use  the  two- 
component  model  given  by  Equations  (19)  and  (20)  to  describe 
the  non-Gaussian  distribution  and  velocity  shear  seen  for  disk 
stars.  Furthermore,  we  assume  that  vR  and  vz  have  uncorrelated 
Gaussian  distributions,  with  zero  mean  and  the  Z-dependent 
intrinsic  dispersion  parameters  listed  in  Table  1.  As  discussed 
in  Section  4,  there  is  no  compelling  evidence  for  a  tilt  in  the 
velocity  ellipsoid  of  blue  disk  stars  in  the  vR-vz  plane,  so  we 
model  the  disk  velocity  ellipsoid  in  cylindrical  coordinates. 


31  They  extrapolated  the  mean  azimuthal  motion  of  color-selected  samples, 
which  is  correlated  with  the  radial-velocity  dispersion,  to  zero  dispersion  and 
obtained  u® ,pec  =  —5.3  km  s-1,  used  here. 


18 


BOND  ET  AL. 


Vol.  716 


-200-100  0  100  200 


ISO" 


metol-poor.  2<D/kpc<7,  dole  <vrod> 


metol-poor.  2<D/kpc<7.  residuol  (D-M)  for  <v„4> 


Figure  16.  Comparison  of  medians  and  dispersions  for  the  measured  and  modeled  radial  velocities  of  20,000  blue  (0.2  <  g  —  r  <  0.4)  halo  stars  (spectroscopic 
[Fe/H]  <  —1.1)  at  distances,  D  =  2  —  7  kpc,  and  b  >  20°.  The  top-left  panel  shows  the  median  measured  radial  velocity  in  each  pixel,  color-coded  according  to  the 
legend  shown  at  the  top  (units  are  km  s-1).  The  top-right  panel  shows  the  difference  between  this  map  and  an  analogous  map  based  on  model-generated  values  of 
radial  velocity,  using  the  same  scale  as  in  the  top-left  panel.  The  bottom-left  panel  shows  the  dispersion  of  measured  radial  velocities,  color-coded  according  to  the 
legend  above  it.  The  bottom-right  panel  shows  the  ratio  of  this  map  and  an  analogous  map  based  on  model-generated  values  of  radial  velocity,  color-coded  according 
to  the  legend  above  it.  When  the  sample  is  divided  into  1  kpc  distance  shells,  the  behavior  is  similar. 


5.2.  Global  Model  Tests 

Our  model  predicts  distributions  of  the  three  measured 
kinematic  quantities,  ura£i,  IM-,  and  /z^,  for  an  arbitrary  control 
volume  defined  by  color,  magnitude,  and  sky  coordinates.  We 
can  test  the  consistency  of  this  model,  m,  with  our  data,  d ,  by 
computing  the  residuals  of  each  kinematic  quantity, 


X  = 


m 


K + 


1/2’ 


(26) 


where  crd  is  the  measurement  error  of  the  data  and  am  is  the 
dispersion  predicted  by  the  model.  For  all  three  kinematic 
quantities,  we  find  x  —  0  to  within  0.05  for  all  three  quantities, 
with  dispersions  of  ~  1.05-1.1.  While  this  result  is  a  necessary 
condition  for  the  model  to  be  acceptable,  it  is  not  sufficient.  A 
stronger  test  of  the  model  was  shown  in  Figure  9  for  disk  stars, 
but  in  what  follows,  we  perform  additional  tests  that  cover  all 
sky  regions  with  available  data. 

5.2.7.  Tests  With  the  Radial-Velocity  Sample 

Figure  16  compares  the  medians  and  dispersions  for  the  mea¬ 
sured  and  modeled  radial  velocities  of  halo  stars.  The  alignment 
of  the  velocity  ellipsoid  with  spherical  Galactocentric  radius 
is  clearly  seen  in  the  bottom-left  panel,  where  the  data  show 
a  dispersion  gradient  moving  away  from  the  Galactic  Center. 


The  increased  dispersion  toward  /  =  180°  was  misinterpreted 
by  Ivezic  et  al.  (2006)  to  be  a  sign  of  substructure.  There  are 
no  large  discrepancies  between  the  measured  and  predicted  be¬ 
havior — the  median  value  of  the  difference  between  observed 
and  modeled  values  is  4.6  km  s-1,  with  a  scatter  of  19  km  s-1 
(see  top-right  panel).  A  similar  scatter  is  obtained  between  two 
model  realizations  with  the  same  number  of  stars  and  measure¬ 
ment  errors.  The  dispersion  ratio,  erf^/cr™^  is  centered  on  1.13, 
with  a  scatter  of  0.2  (see  bottom-right  panel).  For  pairs  of  model 
realizations,  the  ratio  is  always  centered  on  one  to  within  0.02, 
with  a  scatter  of  0.2,  suggesting  that  the  observed  velocity  dis¬ 
persion  is  about  10%  larger  than  predicted  by  our  smooth  model. 
If  the  residuals  were  due  to  halo  rotation,  we  would  expect  a 
spatial  coherence.  Similarly,  for  disk  stars  at  1  <  d  <  2.5  kpc, 
the  median  radial- velocity  residual  is  2.8  km  s-1,  with  a  scatter 
of  6.6  km  s-1  (not  shown). 

5.2.2.  Tests  Based  on  the  Proper-motion  Sample 

The  large  size  of  the  proper-motion  sample  enables  a  much 
higher  spatial  resolution  when  searching  for  structure  in  the 
model  residuals.  We  have  compared  the  observed  and  modeled 
proper-motion  distributions  in  narrow  bins  of  distance,  across 
the  sky,  and  separately  for  disk  and  halo  subsamples.  As  an 
illustration,  Figures  17  and  18  show  the  median  longitudinal  and 
latitudinal  proper  motion  observed  for  disk  stars.  There  is  very 
little  change  in  the  proper-motion  distribution  among  different 
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Figure  17.  Distribution  of  the  median  longitudinal  proper  motion  in  a  Lambert  projection  of  the  North  Galactic  cap  for  high-metallicity  ([Fe/H]  >  —0.9),  blue 
(0.2  <  g  —  r  <  0.4)  stars,  in  several  distance  bins  (top  left:  58,000  stars  with  D  =  0.8-1. 2  kpc;  top  right:  119,000  stars  with  D  =  1. 8-2.2  kpc;  bottom  left:  72,000 
stars  with  D  =  2. 8-3.2  kpc;  bottom  right:  43,000  stars  with  D  =  3. 8-4.2  kpc).  All  maps  are  color-coded  using  the  same  scale,  shown  in  the  middle  (units  are  mas 
yr-1).  Note  that  the  magnitude  of  the  proper  motion  does  not  change  appreciably  as  the  distance  varies  from  ~2  kpc  to  ~4  kpc;  this  is  due  to  a  vertical  gradient  in  the 
rotational  velocity  for  disk  stars  (see  Figure  3). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


distance  bins,  due  to  the  nearly  linear  vertical  rotational- velocity 
gradient.  The  residuals  for  longitudinal  proper  motion  are  shown 
in  Figure  19.  They  provide  weak  evidence  for  either  substructure 
or  a  radial  gradient  that  is  not  modeled,  but  it  is  difficult  to 
distinguish  between  these  two  possibilities  with  these  data  alone. 

A  comparison  of  disk  and  halo  subsamples  in  a  distance  bin 
centered  on  d  =  4  kpc  is  shown  in  Figure  20.  The  largest  data 
versus  model  discrepancy  for  halo  stars,  seen  in  the  bottom-left 
panel  toward  /  ~  0°,  is  also  seen  from  a  different  viewing  angle 
in  the  top-left  panel  of  Figure  8  (R  ~  6  kpc  and  Z  ~  2  kpc).  It 
is  likely  that  this  discrepancy  is  due  to  contamination  of  the  halo 
sample  by  metal-poor  disk  stars.  Figure  21  shows  the  residuals 
for  halo  stars  selected  from  the  8-10  kpc  distance  bin.  The 
residuals  for  both  proper-motion  components  exhibit  similar 
morphology  to  the  systematic  proper-motion  errors  plotted  in  the 
two  left  panels  of  Figure  1 .  In  this  distance  bin,  they  correspond 
to  velocity  errors  of  ~30  km  s_1 ;  as  such,  kinematic  substructure 
in  the  halo  will  be  difficult  to  discern  with  this  sample  at 
>10  kpc.  We  note  that  it  is  tempting  to  associate  the  coherent 
ill  residuals  toward  /  ~  300°  and  l  ~  60°  with  the  Virgo 
overdensity  (see  J08;  An  et  al.  2009).  However,  the  top-left  panel 
in  Figure  1  clearly  shows  systematic  proper-motion  errors  of  the 
required  amplitude  (~1.5  mas  yr-1)  in  the  same  sky  region. 

Schlaufman  et  al.  (2009)  use  SDSS  radial-velocity  measure¬ 
ments  for  metal-poor  turnoff  stars  to  search  for  pieces  of  cold  de¬ 
bris  streams  (which  they  term  Elements  of  Cold  Halo  substruc¬ 


ture,  or  ECHOs).  In  137  lines  of  sight  they  detect  10  ECHOs. 
The  six  northern  ECHOs  from  their  Table  2  (class  I  peak  detec¬ 
tions)  are  shown  as  white  circles  in  the  bottom  two  panels  of 
Figure  21.  It  seems  plausible  that  four  of  these  ECHOs  might 
be  associated  with  the  Monoceros  stream  ( b  ~  30°),  while  one 
of  the  two  remaining  detections  (at  l  =  162? 4  and  l  =  59? 2) 
is  associated  with  the  Grillmair  &  Dionatos  (2006)  stream.  The 
ECHO  at  l  =  100?7  and  b  =  56?8  remains  unassociated  with 
any  known  substructure. 

We  conclude  that  our  model  reproduces  the  first  and  second 
moments  of  the  velocity  distributions  reasonably  well  for  both 
disk  and  halo  stars.  Except  in  the  region  close  to  the  Monoceros 
stream,  the  non-Gaussian  distribution  for  disk  stars  is  also 
well-described.  On  average,  the  model  agrees  with  the  data  to 
within  ~1  mas  yr-1  for  proper  motions,  and  ~10  km  s-1  for 
radial  velocities. 

5.3.  Constraints  on  flsr  from  Large-scale  Halo  Kinematics 

The  proper-motion  distribution  for  halo  stars  toward  the 
Galactic  poles  depends  only  on  the  difference  between  the 
velocity  of  the  local  standard  of  rest  (given  that  the  solar  peculiar 
motion  is  known  to  ~1  km  s-1),  flsr,  and  v $  for  halo  stars.  At 
least  in  principle,  samples  that  extend  over  a  large  sky  area  can 
be  used  to  provide  constraints  on  both  the  halo  rotation  and  flsr- 
Figure  22  compares  the  radial  velocity  and  longitudinal  proper- 
motion  residuals  between  two  models  with  (i^al°,  flsr)  = 
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Figure  18.  Analogous  to  Figure  17,  except  that  the  latitudinal  proper  motion  is  shown. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


(—20,  180)  km  s  \  and  (20,  220)  km  s  1  (for  both  we  fixed 
flsr  —  Fjal°  =  200  km  s_1,  to  make  the  models  agree  with  the 
data  toward  the  North  Galactic  Pole;  see  Figure  5).  In  order  to 
distinguish  these  two  models  observationally,  systematic  errors 
in  the  radial  velocities  must  be  below  10  km  s_1 ,  and  systematic 
errors  in  the  observed  proper  motions  must  be  below  1  mas  yr_1 . 
These  requirements  are  comparable  to  the  systematic  errors  in 
our  data  set,  so  we  can  only  state  that  the  data  are  consistent 
with  i>lsr  ~  200  km  s-1  (or  similarly,  no  halo  rotation)  to 
an  accuracy  of  ~20  km  s_1.  The  proper-motion  measurements 
from  the  Gaia  survey  (Perryman  et  al.  2001)  will  be  sufficiently 
accurate  to  exploit  the  full  potential  of  this  method. 

5.4.  The  Kinematic-Parallax  Relation 

Any  constraint  on  our  best-fit  model  parameters  that  uses 
proper  motions  is  sensitive  to  the  systematic  errors  in  the  dis¬ 
tance  scale  obtained  from  our  photometric-parallax  relationship. 
However,  since  our  model  is  a  good  fit  to  the  radial-velocity  data 
alone  (see  Figure  16),  it  is  possible  to  estimate  distance  errors  by 
minimizing  the  differences  between  the  observed  and  modeled 
proper-motion  distributions.  Such  a  kinematic-parallax  relation, 
derived  from  a  combination  of  radial  velocity  and  proper-motion 
data  sets,  was  proposed  for  the  solar  neighborhood  by  Binney  & 
Tremaine  (1987).  Unlike  in  the  solar  neighborhood,  our  volume 
is  sufficiently  large  that  kinematics  vary  with  position.  Never¬ 
theless,  it  is  conceptually  the  same  method;  a  dipole  (in  our  case 
a  more  complex  angular  function)  is  fitted  to  the  radial  veloc¬ 
ity  and  proper-motion  distributions,  and  the  ratio  of  the  best-fit 
dipole  magnitudes  constrains  the  distance  scale. 


Using  only  low-metallicity  halo  stars,  we  obtain  a  distance- 
scale  error  of  ~5% — our  adopted  absolute  magnitudes  should 
be  ~0. 1  mag  brighter.  This  offset  is  consistent  with  the  expected 
systematic  errors  in  the  calibration  of  the  photometric-parallax 
relation  (see  108).  In  other  words,  the  adopted  distance  scale 
properly  connects  the  radial  and  proper-motion  distributions. 
This  method  provides  a  much  weaker  distance-scale  constraints 
for  disk  stars  than  for  halo  stars,  because  the  vertical  velocity 
gradient  in  the  disk  leads  to  a  degeneracy  between  errors  in 
the  distance  scale  and  errors  in  the  adopted  disk  velocity.  The 
radial- velocity  measurements  do  not  constrain  the  velocity  scale 
because  the  reference  point  depends  on  distance. 

6.  SUMMARY  AND  DISCUSSION 

This  is  the  first  analysis  based  on  SDSS  data  that  simultane¬ 
ously  studies  the  kinematics  of  the  halo  and  disk  populations. 
Past  studies  of  halo  stars  alone  were  performed  by  Sirko  et  al. 
(2004),  Allende  Prieto  et  al.  (2006),  Carollo  et  al.  (2007),  and 
Smith  et  al.  (2009),  while  disk  samples  ranging  from  nearby  M 
stars  to  distant  F/G  stars  have  been  studied  by  Bochanski  et  al. 
(2007b),  West  et  al.  (2008),  and  Fuchs  et  al.  (2009).  Through¬ 
out  this  paper,  we  have  quantified  the  probability  distribution 
function  /?3(/z/,  /Xb,  Fradl u  —  g,  g,  g  —  r,  l,  b),  introduced  in  108, 
that  describes  proper-motion  and  radial-velocity  measurements 
in  the  g  versus  g  —  r  color-magnitude  diagram  as  a  function  of 
position  on  the  sky  and  u  —  g  color.  We  have  developed  a  simple 
empirical  model  with  disk  and  halo  components  that  map  well 
to  populations  detected  in  the  stellar  density  distribution  (J08) 
and  the  metallicity  distribution  (108).  At  distances  accessible 
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metol-rich,  0.8  <  D/kpc  <  1.2,  <resi<Juols>  metol-rich,  1.8  <  D/kpc  <  2.2,  <residuols> 
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Figure  19.  Similar  to  Figure  17,  except  that  the  median  difference  between  the  observed  value  of  longitudinal  proper  motion  and  a  value  predicted  by  the  model 
described  in  the  text  is  shown.  All  maps  are  color-coded  using  the  same  scale,  shown  in  the  middle.  Note  that  the  displayed  scale  is  stretched  by  a  factor  of  2  compared 
to  the  scale  from  Figure  17,  in  order  to  emphasize  discrepancies. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


to  the  HIPPARCOS  survey  (<100  pc),  we  obtain  encouraging 
agreement  with  results  from  Dehnen  &  Binney  (1998),  Dehnen 
(1998),  and  Nordstrom  et  al.  (2004).  The  extension  of  kinematic 
mapping  to  distances  up  to  ~  10  kpc  with  millions  of  stars  repre¬ 
sents  a  significant  observational  advance,  and  delivers  powerful 
new  constraints  on  the  dynamical  structure  of  the  Galaxy.  In 
less  than  two  decades,  the  observational  material  for  such  in 
situ  mapping  has  progressed  from  first  pioneering  studies  based 
on  only  a  few  hundred  objects  (Majewski  1992),  to  over  a  thou¬ 
sand  objects  (Chiba  &  Beers  2000),  to  the  massive  data  set 
discussed  here. 

6.1.  Disk  Kinematics 

The  disk  kinematics  of  the  Milky  Way  are  dominated  by 
rotation  with  a  smooth  vertical  gradient.  Our  analysis  extends 
the  early  measurements  of  this  gradient  (e.g.,  Murray  1986; 
Majewski  1992;  Chiba  &  Beers  2000)  to  vastly  larger  sky 
area  and  to  a  much  larger  distance  range.  The  mean  rotational 
velocity  and  the  three  velocity  dispersions  for  disk  stars  can 
be  modeled  as  simple  functions  of  the  form  a  +  b\Z\c  (see 
Equations  (19)  and  (20),  and  Table  1).  The  rotational- velocity 
distribution  for  the  disk  component  is  non-Gaussian,  and  can 
be  formally  modeled  as  a  sum  of  two  Gaussian  components 
with  fixed  normalization  ratio  for  0.1  <  |Z|/kpc  <  4.  The  fact 
that  the  normalization  ratio  of  these  two  components  does  not 
vary  with  Z  is  at  odds  with  the  standard  disk  decomposition 
into  thin-  and  thick-disk  components  (see  also  Sections  3.4.4 


and  4.2.1  in  108).  Based  on  N-body  simulations  performed  by 
Roskar  et  al.  (2008),  Loebman  et  al.  (2008)  argued  that  the 
absence  of  a  velocity-metallicity  correlation  at  the  thin/thick- 
disk  boundary,  pointed  out  by  108,  may  be  due  to  a  combination 
of  a  strong  vertical  age  gradient  and  the  radial  migration  of  stars 
(see  also  Schonrich  &  Binney  2009).  A  more  detailed  study 
will  be  presented  in  S.  Loebman  et  al.  (2010,  in  preparation).  A 
significant  vertical  age  gradient  for  disk  stars  is  also  supported 
by  an  analysis  of  active  M  dwarfs  presented  in  Bochanski  et  al. 
(2007b).  Such  an  age  gradient,  together  with  the  measured 
velocity  dispersion-age  correlations  for  local  disk  stars  (e.g., 
Nordstrom  et  al.  2004;  Rocha-Pinto  et  al.  2004;  West  et  al. 
2008),  may  be  responsible  for  the  measured  increase  of  velocity 
dispersions  with  distance  from  the  Galactic  plane. 

Close  to  the  plane,  the  proper-motion  data  imply  a  complex 
multi-modal  velocity  distribution  that  is  inconsistent  with  a  de¬ 
scription  based  on  a  simple  Schwarzschild  ellipsoid.  It  is  reas¬ 
suring  that  we  obtained  a  velocity-distribution  morphology  very 
similar  to  that  obtained  by  Dehnen  (1998)  using  statistical  de¬ 
projection  of  the  HIPPARCOS  data,  and  directly  by  Nordstrom 
et  al.  (2004)  using  a  subsample  of  HIPPARCOS  stars  with  full 
three-dimensional  velocity  information.  In  addition,  our  results 
for  the  first  and  second  moments  of  the  velocity  distribution  for 
nearby  M  stars  agree  with  analogous  results  obtained  recently 
by  Fuchs  et  al.  (2009).  The  orientation  of  the  velocity  ellipsoid 
is  strongly  affected  by  multi-modal  structure,  so  one  should  take 
care  in  its  interpretation.  We  discuss  these  issues  in  more  detail 
in  A.  Kowalski  et  al.  (2010,  in  preparation). 
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Figure  20.  Similar  to  Figures  17  and  19,  except  that  the  behavior  of  high-metallicity  (left)  and  low-metallicity  (right)  stars  is  compared  in  a  single  distance  bin 
(3. 5-4. 5  kpc).  The  top  two  panels  show  the  median  longitudinal  proper  motion,  and  the  two  bottom  panels  show  the  median  difference  between  the  observed  and 
model-predicted  values.  An  analogous  figure  for  the  latitudinal  proper  motion  has  similar  characteristics. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


6.2.  Halo  Kinematics 

Our  results  for  the  velocity  distribution  of  halo  stars  are  in 
excellent  agreement  with  Smith  et  al.  (2009).  They  used  proper- 
motion  measurements  based  only  on  SDSS  data,  and  thus  have 
significantly  different,  and  most  likely  much  smaller,  system¬ 
atic  errors  than  the  SDSS-POSS  proper- motion  measurements 
analyzed  here.  The  much  larger  size  of  the  Munn  et  al.  (2004) 
catalog  analyzed  here  allows  us  to  rule  out  the  possibility  that 
the  Smith  et  al.  result  was  biased  by  local  substructure.  The 
close  agreement  of  our  results  for  the  orientation  and  size  of 
the  halo  velocity  ellipsoid  (we  obtained  =  141  km  s-1, 
of  =  75  km  s-1,  and  a £  =  85  km  s_1,  while  their  values 
are  142  km  s_1,  77  km  s_1,  and  81  km  s_1,  respectively)  are 
encouraging  (see  also  Carollo  et  al.  2009).  Their  estimated  er¬ 
rors  of  2  km  s-1  apparently  do  not  include  systematic  effects 
(such  as  errors  in  photometric  parallax;  both  studies  used  the 
same  calibration  from  108) — based  on  our  Monte  Carlo  simu¬ 
lations,  we  believe  that  the  true  errors  cannot  be  smaller  than 
~5  km  s-1.  Additional  independent  evidence  for  the  tilt  of  the 
velocity  ellipsoid  comes  from  the  RAVE  survey  (Siebert  et  al. 
2008).  These  measurements  of  the  velocity  ellipsoid  for  halo 
stars  represent  a  strong  constraint  for  the  shape  of  gravitational 
potential,  as  discussed  by,  e.g.,  Amendt  &  Cuddeford  (1991), 
Kent  &  de  Zeeuw  (1991),  and  Smith  et  al.  (2009). 

We  note  that  Majewski  (1992)  measured  a  retrograde  halo 
rotation  using  stars  observed  toward  the  North  Galactic  Pole  (in 


our  nomenclature,  he  obtained  a  mean  rotational  velocity  v ®  = 
50  zb  16  km  s_1).  A  star-by-star  comparison  of  his  data  and  the 
data  analyzed  here  showed  that  photometric  and  proper-motion 
measurements  agree  within  the  stated  errors.  The  main  reason 
for  different  conclusions  about  halo  rotation  are  the  different 
distance  scales:  his  distances  are  on  average  30%  larger  than  our 
distances,  resulting  in  larger  tangential  velocities.  The  Carollo 
et  al.  (2007)  claim  for  a  large  outer-halo  retrograde  rotation, 
further  refined  by  Carollo  et  al.  (2009)  (in  our  nomenclature,  a 
mean  rotational  velocity  v®  =  80  d=  13  km  s-1),  remains  intact, 
as  their  distance  scale  is  similar  to  ours. 

The  kinematic  measurements  for  halo  stars  presented  here 
should  not  be  extrapolated  beyond  the  sample  distance  limit 
of  10  kpc.  For  example,  using  241  halo  objects,  including 
stars,  globular  clusters,  and  satellite  galaxies,  Battaglia  et  al. 
(2005)  detected  a  continuous  decline  of  the  radial- velocity 
dispersion  beyond  a  Galactocentric  radius  of  ~30  kpc,  from 
about  120  km  s-1  to  50  km  s-1  at  ~120  kpc.  In  addition, 
the  distance  limit  of  our  sample,  together  with  the  decreasing 
sensitivity  of  the  photometric-metallicity  indicator  for  [Fe /H]  < 
—2.5,  prevents  us  from  robustly  testing  the  possible  halo 
dichotomy  discussed  by  Carollo  et  al.  (2007,  2009). 

6.3.  Kinematic  Substructure 

The  model  developed  here  can  be  used  to  search  for  kinematic 
substructure  with  a  low  contrast  level.  For  example,  Schlaufman 
et  al.  (2009)  had  to  generate  a  background  model  when  searching 
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metol-poor.  8  <  D/kpc  <10.  </*l> 


metol-poor.  8  <  D/kpc  <10.  <fit> 
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Figure  21.  Similar  to  Figure  20,  except  that  the  behavior  of  low-metallicity  stars  from  the  8-10  kpc  distance  bin  is  analyzed.  The  top  two  panels  show  the  median 
longitudinal  (left)  and  latitudinal  (right)  proper  motions,  and  the  two  bottom  panels  show  the  median  difference  between  the  observed  and  model-predicted  values.  The 
maps  are  color-coded  according  to  the  legends  in  the  middle  (mas  yr-1;  note  that  the  bottom  scale  has  a  harder  stretch  to  emphasize  structure  in  the  residual  maps). 
The  two  bottom  panels  display  very  similar  morphology  to  systematic  proper-motion  errors  shown  in  the  two  left  panels  in  Figure  1 .  In  the  bottom  panels,  the  white 
symbols  show  the  positions  of  the  six  northern  cold  substructures  identified  by  Schlaufman  et  al.  (2009). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


metol-poor,  2<D/kpc<7.  resid.  (Ml -M2)  for  <vf0d>  metol-poor.  2<D/kpc<7,  doto  <|i,“1  -  2> 


Figure  22.  Comparison  of  radial-velocity  residuals  (left  panel;  analogous  to  the  top-right  panel  in  Figure  16)  and  longitudinal  proper-motion  residuals  (right  panel; 
analogous  to  the  bottom-left  panel  in  Figure  20,  except  for  the  larger  distance  range)  for  two  halo  models  with  (i^al°,  ulsr)=  (—20, 1 80)  km  s_1 ,  and  (20,  220)  km  s_1 . 
Note  that  we  set  i>lsr  —  =  200  km  s_1 .  The  residuals  are  color-coded  according  to  the  legend  above  each  panel  (units  are  km  s-1  for  the  left  panel  and  mas  yr-1 

for  the  right  panel).  In  order  to  distinguish  these  models,  systematic  errors  in  radial  velocity  must  be  below  10  km  s-1,  and  systematic  errors  in  proper  motion  must 
be  below  1  mas  yr-1. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


for  cold  streams  using  radial- velocity  data;  a  similar  study  for 
the  solar  neighborhood  was  performed  by  Klement  et  al.  (2009). 
A  global  model-based  description  is  especially  important  when 


using  large  numbers  of  proper-motion  measurements  to  search 
for  substructure.  While  radial-velocity  data  are  superior  at  large 
distances,  searches  based  on  proper  motions  should  be  better 
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within  a  few  kpc,  due  to  the  high  completeness  and  much 
larger  sample  size.  A  user-friendly  interface  to  our  model  code 
Galfast ,32  which  allows  generation  of  mock  catalogs  in  an 
arbitrary  direction  (or  across  the  entire  sky)  and  to  an  arbitrary 
depth,  will  be  described  elsewhere  (M.  Juric  et  al.  2010,  in 
preparation). 

The  Monoceros  stream  is  clearly  detected  as  a  major  outlier 
from  the  smooth  model  presented  here.  We  have  not  found  any 
other  large  kinematic  substructure  within  10  kpc  that  deviates 
at  a  detectable  level.  Other  deviations  from  the  smooth  model 
predictions  are  likely  due  to  systematic  errors  in  the  proper- 
motion  and  radial-velocity  measurements.  However,  since  the 
main  goal  of  this  paper  was  to  quantify  the  overall  kinematic 
behavior,  the  emphasis  of  our  analysis  was  on  the  first  and 
second  moments  of  the  kinematic  quantities.  It  is  likely  that 
more  sophisticated  statistical  methods,  such  as  those  discussed 
by  Schlaufman  et  al.  (2009)  and  Klement  et  al.  (2009),  will 
be  more  efficient  in  searching  for  substructure.  As  an  example, 
we  have  verified  that  a  moving  group  discovered  by  Majewski 
(1992)  is  reproducible  with  our  data. 

6.4.  Future  Surveys 

The  results  presented  here  will  be  greatly  extended  by  several 
upcoming  large-scale,  deep  optical  surveys,  including  the  Dark 
Energy  Survey  (Flaugher  2008),  Pan-STARRS  (Kaiser  et  al. 
2002),  and  the  Large  Synoptic  Survey  Telescope  (LSST;  Ivezic 
et  al.  2008b).  These  surveys  will  extend  the  faint  limit  of 
this  sample  and  that  of  the  upcoming  Gaia  mission  (Perryman 
et  al.  2001;  Wilkinson  et  al.  2005)  by  4-6  mag.  For  example, 
LSST  will  obtain  proper-motion  measurements  of  comparable 
accuracy  to  those  of  Gaia  at  their  faint  limit,  and  smoothly  extend 
the  error  versus  magnitude  curve  deeper  by  5  mag  (for  details  see 
L.  Eyer  et  al.  2010,  in  preparation).  With  its  w-band  data,  LSST 
will  enable  studies  of  metallicity  and  kinematics  using  the  same 
stars  out  to  a  distance  of  ~100  kpc  (~200  million  F/G  main- 
sequence  stars  brighter  than  g  =  23.5,  for  a  discussion  see  108). 
By  comparison,  the  best  measurements  of  the  outer-halo  radial- 
velocity  dispersion  to  date  are  based  on  from  several  hundred 
(Battaglia  et  al.  2005)  to  several  thousand  (Xue  et  al.  2008) 
objects.  These  upcoming  studies  are  thus  certain  to  provide 
valuable  new  information  about  the  formation  and  evolution  of 
our  Galaxy. 
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APPENDIX  A 

THE  REVISED  SDSS  METALLICITY  SCALE 

Analysis  of  the  metallicity  and  kinematics  of  halo  and  disk 
stars  by  108  utilized  photometric-metallicity  estimates  for  F/G 
stars  with  0.2  <  g  —  r  <  0.6.  Their  mapping  function  from  the 
g  —  r  versus  u  —  g  color-color  diagram  to  metallicity  was  cal¬ 
ibrated  using  stars  with  spectroscopic  metallicities  distributed 
in  SDSS  Data  Release  6.  At  that  time,  high-metallicity  stars 
required  for  the  calibration  of  methods  implemented  in  the 
automated  spectroscopic  pipeline  (SEGUE  Stellar  Parameters 
Pipeline;  Beers  et  al.  2006)  were  not  available.  Between  Data 
Releases  6  and  7,  the  required  data  were  collected  and  the  new 
calibration  resulted  in  the  revised  spectroscopic-metallicity  val¬ 
ues  distributed  with  Data  Release  7  (Lee  et  al.  2008a,  2008b; 
Allende  Prieto  et  al.  2008). 

Here,  we  recalibrate  the  photometric-metallicity  estimator 
using  updated  spectroscopic  metallicities  from  Data  Release  7. 
In  addition,  we  re-derive  the  parts  of  the  108  analysis  that  are 
most  affected  by  this  change  in  the  metallicity  scale. 


A.l.  The  Updated  Photometric -Metallicity  Estimator 

As  shown  in  Figure  A.l,  the  largest  difference  between  the 
SDSS  spectroscopic-metallicity  values  distributed  with  Data 
Releases  6  and  7  is,  as  expected,  at  the  high-metallicity  end. 
In  particular,  the  abrupt  cutoff  in  the  metallicity  distribution  at 
[Fe/H]  ~  —0.5  (see  Figure  9  in  108)  is  no  longer  present,  and 
the  distribution  extends  to  values  as  high  as  [Fe/H]  ~  —  0.2  (the 
distances  for  the  stars  shown  range  from  ~1  kpc  to  ~7  kpc). 

We  proceed  to  re-derive  the  photometric-metallicity  calibra¬ 
tion  using  the  same  selection  criteria  and  the  same  methodology 
as  in  108.  The  new  data  set  admits  a  slightly  simpler  func¬ 
tion — the  double  definition  of  the  v-axis  is  no  longer  required, 
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Figure  A.l.  Summary  of  the  differences  between  SDSS  spectroscopic-metallicity  values  distributed  with  Data  Releases  6  and  7.  The  left  panel  shows  the  median 
difference  between  the  DR7  and  DR6  values  for  0.02  x  0.02  mag2  bins  in  the  g  —  r  vs.  u  —  g  color-color  diagram,  color-coded  according  to  the  legend  shown  in 
the  panel.  The  largest  differences  of  0.2-0.3  dex  are  seen  in  the  top-right  corner,  which  corresponds  to  high  metallicities.  The  right  panel  shows  the  difference  in 
metallicities  as  a  function  of  the  new  DR7  values.  Individual  stars  are  shown  as  small  dots,  and  the  median  values  of  the  difference  are  shown  as  large  circles.  The  two 
dashed  lines  mark  the  ±2cr  envelope  around  the  medians,  where  cr  is  the  rms  scatter  (~0. 1  dex,  due  to  software  updates)  estimated  from  the  interquartile  range.  The 
median  differences  are  larger  than  0.1  dex  only  at  the  high-metallicity  end  ([Fe/H]  >  —0.6). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


and  the  new  expression  is 

[Fe/H]ph  =  A  +  Bx  +  Cy  +  Dxy  +  Ex 2  +  Fy 2  +  Gx2y 

+  Hxy2  +  lx3  +  Jy3,  (Al) 

with  x  =  (u  —  g)  and  y  =  (g  —  r).  The  best-fit  coefficients 
are  (A-J)  =  (-13.13,  14.09,  28.04,  -5.51,  -5.90,  -58.68, 
9.14,  —20.61, 0.0,  58.20).  Note  that  the  coefficient  forx3,  /,  is  to 
zero.  We  removed  this  term  because  it  was  producing  too  much 
curvature  at  the  right  end  (red  u  —  g)  of  the  best-fit  map. 

We  estimate  that  an  upper  limit  on  the  intrinsic  metallicity 
scatter  for  fixed  noiseless  u  —  g  and  g  —  r  colors  (presumably 
due  to  limited  sensitivity  of  broadband  colors  to  metallicity 
variations)  is  about  0.1  dex.  This  value  is  estimated  from  the 
scatter  in  the  difference  between  spectroscopic  and  photometric 
metallicities,  discussed  below.  Unlike  108,  who  simply  adopted 
the  median  metallicity  value  given  by  the  above  expression  for 
each  star,  we  draw  photometric-metallicity  estimates  from  a 
Gaussian  distribution  centered  on  the  best-fit  median  value,  and 
with  a  width  of  0. 1  dex.  The  main  benefit  to  this  is  that  we  avoid 
hard  edges  in  the  photometric-metallicity  distribution  for  stars 
close  to  the  edges  of  the  calibration  region  in  the  g  —  r  versus 
u  —  g  diagram. 

The  performance  of  the  new  map  is  qualitatively  similar 
to  that  of  the  old  map.  The  mean  and  rms  values  for  the 
difference  between  spectroscopic  and  photometric  metallicities 
as  functions  of  the  g  —  r  and  u  —  g  colors  are  shown  in  the  top 
two  panels  of  Figure  A.2.  Typical  systematic  errors  in  the  map 
(i.e.,  the  median  difference  per  pixel)  are  ~0.1  dex  or  smaller, 
and  the  scatter  varies  from  ~0.2  dex  at  the  high-metallicity  end 
to  ~0.3  dex  at  the  low-metallicity  end  (note  that  this  scatter 
includes  a  contribution  from  errors  in  both  spectroscopic  and 
photometric  metallicities). 

The  above  photometric-metallicity  estimator  is  applicable  for 
stars  with  0.2  <  g  —  r  <  0.6  and  — 0.25  +  0.5(w  —  g)  <  g  —  r  < 
0.05  +  0.5 (u  —  g);  that  is,  for  main-sequence  F  and  G  stars  in  the 
calibration  region  of  the  g  —  r  versus  u  —  g  color-color  diagram 
(top  two  panels  of  Figure  A. 2).  For  stars  with  spectroscopic 
metallicity  [Fe/H]  >  —2.2,  the  distribution  of  the  difference 
between  spectroscopic  and  photometric  metallicities  is  well- 
described  by  a  Gaussian  with  a  width  of  0.26  dex  (see  the 
bottom-right  panel  in  Figure  A.2). 


It  should  be  noted  that  the  performance  of  the  photometric- 
metallicity  estimator  deteriorates  at  the  low-metallicity  end  be¬ 
cause  the  u  —  g  color  becomes  insensitive  to  decreases  in  metal¬ 
licity.  As  shown  in  the  bottom-left  panel  of  Figure  A.2,  the 
photometric  metallicity  saturates  at  [Fe/H]  ~  —2  for  smaller 
values  of  spectroscopic  metallicity.  Even  at  [Fe/H]  =  —2, 
the  true  metallicity  is  overestimated  by  0. 2-0.3  dex,  and  by 
[Fe/H]  =  —  3  this  bias  is  as  large  as  1  dex  (the  photometric- 
metallicity  values  never  become  significantly  lower  than 
[Fe/H]  =  —  2).  This  shortcoming  could  be  partially  alleviated 
by  employing  more  accurate  u- band  photometry  (say,  with  er¬ 
rors  of  0.01  mag  instead  of  0.03  mag,  as  used  here),  but  probably 
not  for  metallicities  lower  than  [Fe/H]  =  —2.5.  Fortunately,  the 
low-metallicity  inner-halo  stars  within  SDSS  reach  have  a  me¬ 
dian  metallicity  of  [Fe/H]  ~  — 1.5  (108),  so  for  the  majority  of 
stars  the  photometric  metallicities  are  robust.  Another  important 
note  is  that,  despite  the  improvement  at  the  high-metallicity  end, 
the  calibration  range  only  extends  to  [Fe/H]  ~  —0.2.  Our  cali¬ 
bration  sample  did  not  include  young,  metal-rich  main-sequence 
stars  (w  —  g  >  1.1  and  g  —  r  <  0.3).  For  this  reason,  our  polyno¬ 
mial  model  underestimates  true  metallicities  by  0. 2-0.3  dex  at 
[Fe/H]  =  0,  and  probably  more  for  [Fe/H]  >  0.  Any  result  re¬ 
lying  on  higher  metallicities  should  be  interpreted  with  caution 
(especially  at  low  Galactic  latitudes  where  the  uncertain  in¬ 
terstellar  medium  extinction  may  strongly  affect  the  estimated 
metallicities).  We  plan  to  extend  our  calibration  further  into 
the  metal-rich  domain  by  employing  data  from  the  SDSS-III 
SEGUE-2  survey  (Yanny  et  al.  2009),  and  further  refinements 
of  the  SEGUE  Stellar  Parameters  Pipeline,  which  are  now  un¬ 
derway. 

A.2.  Tomography  II Reloaded 

108  pointed  out  several  aspects  of  their  analysis  that  may 
have  been  affected  by  the  metallicity  “compression”  at  the 
high-metallicity  end  in  DR6.  We  repeated  their  full  analysis 
and  report  here  on  those  aspects  where  differences  warrant 
discussion. 

The  “hard”  upper  limit  on  photometric-metallicity  estimates 
at  the  high-metallicity  end  ([Fe/H]  ~  —0.5)  with  the  DR6 
calibration  is  best  seen  in  the  bottom-left  panel  in  Figure  9 
from  108.  We  reproduce  that  map  of  the  conditional  metallicity 
distribution  in  the  top-left  panel  of  Figure  A. 3.  As  expected,  the 
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Figure  A.2.  Color-coded  map  in  the  top-left  panel  shows  the  median  difference  between  spectroscopic  and  revised  photometric  metallicities  for  ~50,000  stars  from 
SDSS  Data  Release  7.  The  median  value  is  zero,  and  the  rms  scatter  is  0.07  dex.  The  contours  show  the  distribution  of  stars  with  r  <  20  and  at  high  Galactic  latitudes. 
The  top-right  panel  shows  the  rms  scatter  of  the  difference  between  spectroscopic  and  photometric  metallicities  in  each  pixel.  The  top  two  panels  are  analogous  to  the 
bottom  two  panels  in  Figure  2  from  108.  The  bottom-left  panel  shows  the  photometric  metallicity  as  a  function  of  the  spectroscopic  metallicity.  Individual  stars  are 
shown  by  small  dots,  and  the  median  values  of  the  difference  are  shown  by  large  circles.  The  distribution  of  stars  is  shown  as  linearly  spaced  contours.  Note  that  the 
photometric  metallicity  saturates  at  [Fe/H]  ~  —  2  at  the  low-metallicity  end.  The  histogram  in  the  bottom-right  panel  shows  the  distribution  of  the  difference  between 
spectroscopic  and  photometric  metallicities  for  stars  with  spectroscopic  metallicity  [Fe/H]  >  —2.2.  A  best-fit  Gaussian  centered  on  zero  and  with  a  width  of  0.26  dex 
is  shown  by  the  dashed  line. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


metallicity  distribution  of  disk  stars  within  2  kpc  of  the  Galactic 
plane  now  extends  to  [Fe/H]  ~  0. 

In  the  new  calibration,  the  parameters  describing  the  variation 
of  the  median  metallicity  for  disk  stars  as  a  function  of  the 
distance  from  the  Galactic  plane, 

/X/)(Z)  =  [1  oo  +  exp(— |Z|///m)  dex,  (A2) 

are  also  changed.  The  updated  values  are  H ^  =  0.5  kpc,  = 
—0.82,  and  =  0.55  (the  old  values  were  H ^  =  1.0  kpc, 
/Too  =  —0.78,  and  A^  =  0.35).  The  best-fit  values  of 
and  Am  are  accurate  to  ~0.05  dex.  Values  of  H M  in  the  range 
350-700  pc  are  consistent  with  the  data — the  decrease  in  H M 
is  required  by  the  local  constraint,  /7£>(Z)  =  —0.2  (Nordstrom 
et  al.  2004;  Allende  Prieto  et  al.  2004). 

An  interesting  result  from  108  was  the  detection  of  disk  stars 
at  a  distance  from  the  Galactic  plane  as  large  as  ~6  kpc  (see 
their  Figure  10).  A  peak  at  [Fe/H]  =  —0.5  in  the  metallicity 
distribution  of  stars  at  those  distances  was  another  manifestation 
of  the  metallicity  “compression.”  As  demonstrated  in  the  top- 
right  panel  of  Figure  A.3,  this  peak  is  not  present  when  using  the 
revised  calibration.  However,  there  is  still  statistical  evidence 
that  disk  stars  exist  at  these  large  distances  from  the  plane: 
about  5%  of  stars  in  the  5  <  Z  <  7  kpc  bin  are  consistent  with 
disk  stars,  in  agreement  with  extrapolation  of  the  exponential 
profile  derived  from  stellar  counts.  Lee  &  Beers  (2009)  and 
Carollo  et  al.  (2009)  have  also  commented  on  signatures  of 


thin-disk-like  chemistry  and  kinematics  for  a  small  fraction  of 
stars  several  kpc  above  the  Galactic  plane. 

Perhaps  the  most  intriguing  result  of  108  was  the  non¬ 
detection  of  a  correlation  between  rotational  velocity  and 
metallicity  for  disk  stars  at  Z  ~  1  kpc.  At  such  distances 
from  the  Galactic  plane,  the  counts  of  thin-  and  thick-disk 
stars  are  expected  to  be  similar.  Since  traditionally  the  thick- 
disk  component  is  associated  with  a  larger  velocity  lag  and 
lower  metallicities,  a  fairly  strong  and  detectable  correlation 
was  expected  (see  108  for  details).  The  two  bottom  panels  in 
Figure  A.3  demonstrate  that  such  a  correlation  is  still  undetected, 
although  the  photometric-metallicity  range  now  extends  to 
higher  values  (up  to  [Fe/H]  ~  —0.2). 

The  higher  metallicity  values  obtained  with  the  recalibrated 
relation  have  a  quantitative  effect  on  the  best-fit  metallicity  dis¬ 
tributions  shown  in  Figure  7  of  108.  Using  the  same  method¬ 
ology,  we  reproduce  the  metallicity  distributions  with  the  new 
calibration  in  Figure  A.4.  108  modeled  the  non-Gaussian  disk 
metallicity  distribution  using  a  sum  of  two  Gaussians  with  a 
fixed  amplitude  ratio  (1.7:1),  fixed  difference  of  the  mean  val¬ 
ues  (0. 14  dex),  and  fixed  widths  (0.21  dex  and  0.11  dex),  which 
“slides”  as  a  function  of  Z,  according  to  Equation  (A2).  We  find 
that  the  only  significant  change  is  an  increase  to  the  width  of 
the  second  Gaussian  to  0.21  dex,  which  accounts  for  the  exten¬ 
sion  of  the  metallicity  distribution  to  higher  values.  Only  minor 
changes  are  required  for  the  best-fit  halo  metallicity  distribution 
(see  Table  3  in  108):  the  median  halo  metallicity  is  now  —1.46 
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NCP.  0.8  <  Z/kpc  <  1-2  NCP.  1.2  <  2/kpc  <  1.4 


Figure  A.3.  Conditional  metallicity  distribution  as  a  function  of  distance  from  the  Galactic  plane  (top  left),  analogous  to  the  bottom-left  panel  of  Figure  9  from  108. 
Note  that  the  new  photometric  metallicities  include  [Fe/H]  >  —0.4.  The  dashed  line,  which  shows  the  median  disk  metallicity,  is  also  revised  (see  the  text).  The 
top-right  panel  is  analogous  to  Figure  10  from  108,  and  shows  the  metallicity  distribution  for  stars  with  5  <  Z  <  7  kpc,  where  Z  is  the  distance  from  the  Galactic  plane. 
Note  that  the  photometric-metallicity  artifact  at  [Fe/H]  =  —0.5  discussed  by  108  is  no  longer  present.  However,  there  is  still  evidence  that  disk  stars  exist  at  such  large 
distances  from  the  plane.  The  bottom  two  panels  show  the  heliocentric  rotational  velocity  for  disk  stars  in  two  thin  Z  slices,  and  are  analogous  to  the  bottom-right 
panel  in  Figure  16  from  108.  Note  that  the  correlation  between  velocity  and  metallicity  is  still  absent. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Figure  A.4.  Analogous  to  Figure  7  from  108.  The  symbols  with  error  bars  show  the  metallicity  distribution  for  stars  with  0.2  <  g  —  r  <  0.4,  7  <  R  <  9  kpc, 
and  distances  from  the  Galactic  plane  as  marked,  where  R  is  the  Galactocentric  cylindrical  radius.  The  behavior  is  qualitatively  similar  to  that  seen  in  108.  The  only 
significant  quantitative  difference  is  in  the  model  for  the  metallicity  distribution  of  disk  stars  (see  the  text). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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in  the  first  three  bins,  and  —1.56  in  the  most  distant  Z  bin,  and 
its  width  changed  from  0.32  dex  to  0.36  dex  in  the  last  bin.  We 
note  somewhat  less  scatter  of  the  data  points  around  the  best-fit 
functions  with  the  recalibrated  data  set.  To  summarize,  the  re¬ 
vised  best-fit  parameters  that  describe  halo  and  disk  metallicity 
distributions  are  as  follows. 

1.  The  halo  metallicity  distribution  is  spatially  invariant  and 
well-described  by  a  Gaussian  distribution  centered  on 
[Fe/H]  =  —1.46,  with  a  width  oH  —  0.30  dex  (not  includ¬ 
ing  measurement  errors).  For  \Z\  <  10  kpc,  an  upper  limit 
on  the  halo  radial-metallicity  gradient  is  0.005  dex  kpc-1. 

2.  The  disk  metallicity  distribution  varies  with  Z  such  that 
its  shape  remains  fixed,  while  its  median,  fiD,  varies  as 
given  by  Equation  (A2)  (with  best-fit  parameter  values 

=  0.5  kpc,  /Zoo  =  —0.82,  and  AM  =  0.55).  The  shape 
of  the  disk  metallicity  distribution  can  be  modeled  as 

pD(x  =  [Fe/H] | Z)  =  0.63  G[x\/jl  =  a(Z\  o  =  0.2] 

+  0.37  G\x |/z  =  a(Z)  +  0.14,  a  =  0.2],  (A3) 

where  the  position  a  and  the  median  \jld  are  related  via 
a(Z)  —  fiD{Z)  —  0.067  (unless  measurement  errors  are 
very  large). 

We  point  out  that  the  asymmetry  of  the  metallicity  distribution 
for  disk  stars  is  now  less  pronounced  (as  implied  by  the  same 
widths  of  the  two  best-fit  Gaussian  components).  Nevertheless, 
due  to  our  large  sample  size,  the  non-Gaussianity  is  detected 
at  high  significance.  A  remaining  uncertainty  is  the  error 
distribution  of  the  photometric  metallicities,  which  itself  could 
account  for  such  a  deviation  from  Gaussianity.  However,  to  the 
extent  possible  using  a  highly  incomplete  spectroscopic  sample 
(cf.  the  bottom-right  panel  in  Figure  A. 2  and  discussion  in  108), 
we  are  unable  to  quantitatively  explain  the  observed  deviation 
from  Gaussianity  as  an  artifact  of  the  photometric-metallicity 
method. 

APPENDIX  B 
SQL  QUERY  EXAMPEE 

The  following  SQL  query  was  used  to  select  and  down¬ 
load  data  for  all  SDSS  stars  with  spectroscopic  and  proper- 
motion  measurements  (see  http://casjobs.sdss.org/CasJobs). 
SELECT 

round (p . ra,6)  as  ra,  round (p . dec , 6)  as  dec, 
p.run,  p. rerun,  round (p. extinct ion_r, 3)  as  rExt , 
round (p .psfMag_u,2)  as  upsf ,  —  comments  are 
preceded  by  — 

round (p .psfMag_g,2)  as  gpsf ,  —  rounding  up 

round  (p .psfMag_r ,2)  as  rpsf, 

round(p.psfMag_i,2)  as  ipsf , 

round (p .psfMag_z,2)  as  zpsf , 

round (p .psfMagErr_u,2)  as  uErr, 

round (p .psfMagErr_g,2)  as  gErr, 

round (p .psfMagErr_r ,2)  as  rErr, 

round (p .psfMagErr_i ,2)  as  iErr, 

round (p .psfMagErr_z ,2)  as  zErr, 

round(s .pmL,2)  as  pmL,  round(s .pmB,2)  as  pmB, 

pmRaErr , 

t.specObjID,  t. plate,  t.mjd,  t.fiberid,  t.feha, 
t . f ehaerr , 

t.logga,  t.loggaerr,  t.elodierv,  t . elodierverr 
INTO  mydb.pmSpec 


FROM  star  p,  propermotions  s,  sppParams  t, 

specobjall  q 

WHERE 

p. objID  =  s.objID  and  s. match  =1  —  must  have 
proper  motion  and  t.specobjid  =  q.specobjid  and 

q. bestobjid  =  p.objid  and  s.sigra  <  350  and 
s.sigdec  <  350  —  quality  cut  on  pm  and  (p. flags 
&  ;  4295229440 ;  )=  0  —  see  text  for  flag  list  and 
p.psfMag_r  >  14.5  —  avoid  saturation  and 
p.psfMag_r  <  20  —  practical  faint  limit  for  pm 

—  the  end  of  query 

The  following  SQL  query  was  used  to  select  and  down¬ 
load  data  for  all  spectroscopically  confirmed  quasars  with 
proper-motion  measurements  and  redshifts  in  the  range  0.5-2. 5. 
SELECT 

round (p . ra,6)  as  ra,  round (p . dec ,6)  as  dec, 

p. run,  p. rerun,  round (p . extinct ion_r , 3)  as  rExt, 
round (p .psfMag_u,2)  as  upsf, 
round(p.psfMag_g,2)  as  gpsf, 

round (p .psfMag_r ,2)  as  rpsf, 

round(p.psfMag_i,2)  as  ipsf, 

round (p .psfMagErr_u,2)  as  uErr, 

round (p .psfMagErr_g,2)  as  gErr, 

round (p .psfMagErr_r ,2)  as  rErr, 

round (p .psfMagErr_i ,2)  as  iErr, 

round (s .pmL,2)  as  pmL,  round(s .pmB,2)  as  pmB, 

pmRaErr , 

q. specObjID,  q. plate,  q.mjd,  q.fiberlD, 

q.z,  q.zErr,  q.zConf,  q.zWarning,  q.specClass 
INTO  mydb.pmQSO 

FROM  star  p,  propermotions  s,  specobjall  q 
WHERE 

p. objID  =  s.objID  and  s. match  =  1  and 

q. bestobjid  =  p.objid  and  s.sigra  <  350  and 
s.sigdec  <  350  —  per  s.s  and  (p. flags  & 

; 4295229440 ; )  =  0  and  p.psfMag_r  >  14.5  and 
p.psfMag_r  <  20  and  q.z  >  0.5  —  redshift 
limits  and  q.z  <  2.5 

—  the  end  of  query 
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